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ACTIVE  AND  PASSIVE  REMOTE  SENSING  OF  ICE 
FINAL  REPORT 

This  annual  report  covers  research  under  the  sponsorship  of  the  ONR  contract 
N00014-89-J— 1107  from  October  1,  1988  to  September  30,  1992.  We  have  published  44 
journal  and  conference  papers  and  5  student  theses  during  this  period. 

The  theoretical  approach  that  has  been  developed  to  interpret  the  polarimetric 
active  measurements  of  saline  ice  is  a  random  medium  model  using  the  radiative  trsmsfer 
theory.  The  ice  layer  is  described  as  a  host  ice  medium  embedded  with  rtmdomly  dis¬ 
tributed  inhomogeneities,  and  the  underlying  sea  water  is  considered  to  be  a  homogeneous 
halfspace.  Multiple  scattering  effects  are  accounted  for  by  solving  the  radiative  trsmsfer 
equations  numerically.  The  effects  of  random  roughness  at  the  air  -  ice,  and  ice  -  water 
interfaces  are  accounted  for  by  modifying  the  boundEiry  conditions  in  the  radiative  trans¬ 
fer  equations.  Analysis  of  the  model  for  reconstruction  of  sea  ice  parameters  is  made.  An 
optimization  approach  is  used  for  inversion.  The  discrepzmcy  between  the  data  and  the 
results  of  the  forward  model  is  minimized  by  changing  the  inversion  parzuneters  according 
to  a  nonlinear  programming  scheme.  Reconstruction  of  correlation  lengths  in  the  horizon¬ 
tal  and  vertical  dimensions  heis  been  accomplished  using  the  polarimetric  backscattering 
coefficients  at  different  angles  of  incidence  as  input  data.  Effects  of  data  diversity  and 
noise  on  the  reconstruction  of  the  physical  parameters  of  sea  ice  from  the  backscattering 
coefficients  are  being  investigated. 

Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the 
remote  sensing  of  geophysiczil  media.  It  has  been  shown  that  the  third  Stokes  parameter 
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U  of  the  thermal  emission  may  become  large  for  azimuthally  asymmetric  fields  of  abser- 
vation.  In  one  experiment,  values  of  U  of  as  high  as  40  K  were  meeisured  fi:om  a  periodic 
triangular  soil  surface  at  10  GHz.  In  order  to  investigate  the  potential  applicability  of 
passive  polarimetry  to  the  remote  sensing  of  oceain  surface,  a  numerical  study  of  the  po- 
larimetric  thermal  emission  from  randomly  rough  oceem  surface  W2is  performed.  A  Monte 
Carlo  technique  utilizing  an  exact  method  for  calculating  thermal  emission  was  chosen  for 
the  study  to  avoid  any  of  the  limitations  of  the  commonly  used  approximate  methods  in 
rough  surface  scattering. 

In  this  Monte  Carlo  technique,  a  set  of  finite  rough  surface  profiles  in  two  di¬ 
mensions  with  desired  statistics  was  generated.  Each  finite  surface  was  then  extended 
periodically  to  create  a  set  of  infinite  rough  surfaces.  The  polarimetric  thermal  emission 
from  each  surface  of  the  set  was  then  calculated  using  both  the  extended  boundary  con¬ 
dition  method  and  the  method  of  moments.  Finally,  the  results  from  the  set  were  then 
averaged  to  obtain  the  estimate  of  the  polarimetric  brightness  temperatures  for  the  given 
surface  statistics.  The  surface  statistics  chosen  were  intended  to  model  a  wind  perturbed 
ocean  surface  in  the  X  to  Xu  band  microwave  region.  For  this  purpose,  a  power  law 
surface  spectrum  corresponding  to  the  capillary  wave  portion  of  the  ocean  spectrum  was 
chosen  with  rms  surface  heights  corresponding  to  wind  speeds  of  up  to  10  m/s.  The  effects 
of  varying  the  azimuthal  and  polar  looking  angles,  rms  surface  height,  ocean  permittivity, 
exponent  of  the  power  law  spectrum,  and  the  necessary  high  and  low  frequency  cutoffs  of 
the  spectrum  were  investigated  in  this  numerical  study. 

The  results  of  the  study  indicate  that  the  U  parameter  is  sensitive  to  the  azimuthal 
angle  between  the  surface  periodicity  and  the  looking  angle  and  to  the  rms  height  of 
the  surface,  and  that  the  U  parameter  is  fairly  insensitive  to  variations  in  polar  angle, 
permittivity,  surface  power  law  spectrum,  and  surface  spectrum  high  frequency  cutoff. 
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These  properties  give  further  strength  to  the  idea  of  using  the  U  parameter  to  detect  wind 
direction  over  the  ocean. 

We  studied  an  anisotropic  layer  model  with  ellipsoidal  scatterers  for  applications  to 
polarimetric  remote  sensing  of  geophysical  media  at  microwave  frequencies.  The  scattering 
configuration  includes  an  isotropic  layer  covering  an  anisotropic  layer  above  a  homogeneous 
half  space.  The  isotropic  layer  consists  of  remdomly  oriented  spheroids.  The  anisotropic 
layer  contains  ellipsoidal  scatterers  with  preferential  verticeJ  alignment  and  random  az¬ 
imuthal  orientations.  Effective  permittivities  of  the  scattering  media  are  calculated  with 
the  strong  permittivity  theory  extended  to  account  for  the  non-spherical  shapes  and  the 
scatterer  orientation  distributions.  Based  on  the  aneJytic  wave  theory,  the  dyadic  Green’s 
function  for  layer  media  is  used  to  derive  polarimetric  backscattering  coefficients  under  the 
distorted  Born  approximation.  The  ellipsoidal  shape  of  the  scatterers  gives  rise  to  non¬ 
zero  cross-polarized  returns  from  the  untilted  anisotropic  medium  even  in  the  first  order 
approximation.  The  effect  of  rough  interface  is  estimated  by  incoherent  addition  method. 
Theoretical  results  and  experimental  data  are  matched  at  9  GHz  for  thick  first-year  sea  ice 
with  bsire  surface  and  with  snow  cover  at  Point  Barrow,  Alaska.  The  model  is  then  used 
to  study  the  sensitivity  of  the  scattering  coefficients  with  respect  to  correlation  lengths 
representing  the  geometry  of  brine  inclusions  in  the  model  for  sea  ice.  Fin^diy,  polarimet¬ 
ric  signatures  of  bare  and  snow-covered  sea  ice  are  also  simulated  based  on  the  model  to 
investigate  effects  of  different  scattering  mechanisms. 


A  model  to  calculate  the  effective  permittivity  of  saline  ice  under  thermal  variation 
is  presented  in  this  paper.  The  model  includes  multi-phzise  inhomogeneities  with  multiple 
species  characterized  by  orientation,  size,  and  shape  distributions.  The  model  is  then 
used  to  derive  the  effective  permittivity  as  a  function  of  temperature  under  the  strong 
fluctuation  theory  which  is  extended  to  account  for  the  complexity.  The  results  calculated 
from  the  model  are  compared  with  experimental  data  at  4.8  GHz  for  saline  ice  grown  at  the 
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us  Army  gold  Regions  Research  and  Engineering  Laboratory  (CRREL).  The  comparison 
between  measured  and  calculated  complex  permittivities  is  good  for  the  imaginziry  part 
and  the  difference  is  within  10%  for  the  real  part. 

Fully  polarimetric  scattering  of  electromagnetic  waves  &om  snow  and  ice  is  stud¬ 
ied  with  a  multi-layered  random  medium  model  and  applied  to  interpret  experimental 
data  obtained  under  laboratory  controlled  conditions  such  as  CRRELEX.  The  snow  layer 
is  modeled  as  an  isotropic  random  medium.  The  sea  ice  is  described  as  an  anisotropic 
random  medium  due  to  the  nonspherical  shape  of  brine  inclusions.  The  underlying  sea 
water  is  considered  eis  a  homogeneous  half-space.  The  random  media  in  both  layers  are 
characterized  by  three-dimensional  correlation  functions  with  variances  and  correlation 
lengths  corresponding  to  the  fluctuation  strengths  and  the  physical  geometries  of  the  inho¬ 
mogeneities,  respectively.  The  strong  fluctuation  theory  is  used  to  calculate  the  effective 
permittivities  of  the  random  media.  The  distorted  Born  approximation  is  then  employed 
to  obtain  the  covariance  matrix  which  represents  the  fully  polarimetric  scattering  prop¬ 
erties  of  the  snow-ice  media.  It  has  been  shown  that  the  polarimetric  covariance  matrix 
contains  more  information  than  the  conventional  scattering  coefficients  on  the  remotely 
sensed  media. 

In  saline  ice  sheets  under  quiescent  condition,  the  background  ice  grows  in  colum¬ 
nar  form  and  saline  water  is  trapped  between  ice  platelets  in  the  form  of  brine  inclusions 
which  are  usually  ellipsoidal.  The  ice  tends  to  grow  vertically  downward  rendering  the 
ellipsoidal  inclusions  aligned  preferably  in  the  vertical  direction  and  the  crystallographic 
C  axes  parallel  to  the  horizontal  plane.  In  this  case,  the  C  axes  are,  however,  random  in 
azimuthal  direction.  The  strong  fluctuation  theory  is  extended  to  account  for  vertically 
aligned  ellipsoidal  brine  inclusions  with  C  axes  randomly  oriented  in  the  horizontzd  direc¬ 
tion.  The  brine  inclusions  are  described  by  three-dimensional  local  correlation  functions. 
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The  configuration  average  over  the  azimuthal  orientation  angles  is  carried  out  in  the  pro¬ 
cess  of  deriving  the  global  correlation  tensor.  The  distorted  Born  approximation  is  applied 
to  obtain  the  covariance  matrix  for  the  multi-layered  snow-ice  configuration.  The  theo¬ 
retical  results  show  non-zero  cross-polarized  returns  under  the  first-order  distorted  Born 
approximation.  We  have  also  compared  the  results  with  experimental  data  obtained  by 
the  US  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL). 

In  sea  ice,  the  scatterers  can  have  various  shapes,  sizes,  and  permittivities.  We 
have  also  investigated  the  modeling  of  radar  backscatter  Rom  random  media  with  multiple 
scatterer  species.  We  consider  each  type  of  scatterers  as  a  species  which  can  take  on  a  shape, 
size,  and  complex  permittiAdty  different  from  other  species.  The  multiple  species  in  the 
random  medium  are  considered  as  randomly  oriented  ellipsoids  and  described  by  multiple 
three-dimensional  ellipsoidal  local  correlation  functions.  The  Vciriances  and  correlation 
lengths  of  the  correlation  functions  characterize  the  fluctuation  strengths  and  the  physical 
geometries  of  each  species  of  scatterers.  The  effective  permittivity  of  the  random  medium 
is  derived  under  the  strong  fluctuation  theory  and  the  polarimetric  scattering  coefficients 
are  calculated  for  the  layer  configuration  with  the  distorted  Born  approximation.  Due 
to  the  non-spherical  shape  and  the  random  orientation  of  the  scatterers,  the  correlation 
coefficient  between  the  HH  and  VV  returns  has  a  magnitude  different  from  unity  and  a 
small  phase  angle.  The  scattering  coefficients  are  also  used  to  calculate  the  Mueller  matrix 
for  synthesis  of  polarization  signatures.  The  co-polarized  signature  of  the  random  medium 
heis  a  rather  straight  distortion  track  and  a  recognizable  pedestal. 

Accurate  calibration  of  polarimetric  radar  systems  is  essentied  for  the  polarimetric 
remote  sensing  of  earth  terrain.  Polarimetric  calibration  algorithms  using  in-scene  re¬ 
flectors  and  distributed  targets  eire  studied.  The  transmitting  and  receiving  ports  of  the 
polarimetric  radar  are  modeled  by  two  unknown  polarization  transfer  matrices.  These 
unknown  matrices  are  determined  using  the  the  measured  scattering  matrices  from  the 
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calibration  targets.  For  the  czise  of  polzurimetric  calibration  using  three  in-scene  reflec¬ 
tors,  a  Polarization-basis  Transformation  (PT)  technique  is  introduced  to  find  out  a  new 
transmitting  and  receiving  polarization  basis  under  which  the  scattering  matrices  of  the 
calibration  targets  will  fall  into  one  of  six  simpler  sets.  The  calibration  solution  can  then 
be  solved  easily  in  the  new  polarization  basis  and  converted  to  obtain  the  solution  in  the 
original  polarization  basis.  The  uniqueness  of  polarimetric  calibration  using  three  targets 
is  addressed  for  all  possible  target  combinations.  The  PT  technique  can  also  be  applied 
to  the  polarimetric  calibration  using  a  combination  of  arbitrary  in-scene  reflectors  and 
distributed  targets.  The  effect  of  misalignment  of  calibration  targets  and  the  sensitivity 
of  polarimetric  calibration  algorithms  to  the  noise  are  illustrated  by  investigating  several 
sets  of  calibration  targets. 

We  have  studied  the  SAR  image  classification  by  using  the  neural  network  meth¬ 
ods.  Supervised  methods,  including  both  conventional  Maximum  Likelihood  (ML)  and 
more  recent  multi-layer  perceptron  neural  network  classifiers  have  yielded  higher  accuracy 
than  unsupervised  techniques,  but  suffer  from  the  need  for  human  interaction  to  prede¬ 
termine  classes  and  training  regions.  In  contrast,  unsupervised  methods  determine  classes 
automatically,  but  generally  show  limited  ability  to  accurately  divide  terrain  into  natural 
classes.  We  introduced  a  new  terrsdn  classification  technique  to  discriminate  sea  ice  sig¬ 
natures  in  polarimetric  SAR  images  by  utilizing  unsupervised  neural  networks  to  provide 
automatic  classification,  but  employing  an  iterative  algorithm  which  overcomes  the  poor 
accuracy  of  other  unsupervised  techniques. 

Several  types  of  unsupervised  neural  networks  are  first  applied  to  the  classification 
of  SAR  images,  and  the  results  are  compared  with  those  of  more  conventional  unsupervised 
methods.  Neural  network  approaches  include  Adaptive  Resonance  theory  (ART),  Learning 
Vector  Quantization  (LVQ),  and  Kohonen’s  self-organizing  feature  map.  Conventional 
classifiers  utilized  are  the  migrating  means  clustering  algorithm  and  the  K-means  clustering 
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method.  Preprocessing  is  performed  with  the  SAR  images  to  reduce  speckle  noise  and 
stabilize  the  training  process  for  the  neural  networks.  Results  after  preprocessing  show  that 
LVQ  and  Kohonen’s  self-organizing  feature  map  outperform  the  conventionai  unsupervised 
classifiers,  but  are  still  inferior  to  supervised  methods. 

To  overcome  this  poor  accuracy,  an  iterative  algorithm  is  constructed  where  the 
SAR  image  is  reclassified  using  a  Maximum  Likelihood  (ML)  classifier.  Training  of  the  ML 
classifier  is  performed  using  a  training  data  set  first  classified  by  the  above  unsupervised 
method,  thus,  requiring  no  human  intervention,  and  preserving  the  unsupervised  nature 
of  the  overall  classification  scheme.  The  process  is  then  repeated  iteratively,  training  a 
second  ML  classifier  using  data  classified  by  the  first.  It  is  shown  that  this  algorithm 
converges  rapidly,  and  significantly  improves  claissification  accuracy.  Performance  after 
convergence  is  seen  to  be  comparable  to  that  obtained  with  a  supervised  ML  classifier, 
while  maintaining  the  advantages  of  an  unsupervised  technique. 

The  unsupervised  and  iterative  techniques  developed  have  been  applied  to  the 
polarimetric  SAR  images  of  Beaufort  sea  ice  acquired  by  the  G-,  L-,  and  P-band  SAR 
instruments  of  Jet  Propulsion  Laboratory.  The  results  obtained  with  the  new  algorithms 
are  compared  with  the  resiilts  obtedned  with  other  techniques  by  classifying  terrain  features 
in  polarimetric  SAR  images. 

A  multivariate  K-distribution  has  been  developed  to  model  the  statistics  of  fully 
polarimetric  radar  data  from  earth  terrain  Avith  polarizations  HH,  HV,  VH,  and  VV.  In 
this  approach,  correlated  polarizations  of  radar  signals,  as  characterized  by  a  covariance 
matrix,  are  treated  as  the  sum  of  N  n-dimensional  random  vectors;  N  obeys  the  negative 
binomial  distribution  with  a  parameter  a  and  mean  N.  Subsequently,  an  n-dimensional 
K-distribution,  with  either  zero  or  nonzero  mean,  is  developed  in  the  limit  of  infinite  N  or 
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illuminated  area.  The  probability  density  function  (PDF)  of  the  K-distributed  vector  nor¬ 
malized  by  its  Euclidean  norm  is  independent  of  the  parameter  a  and  is  the  same  as  that 
derived  from  a  zero-mean  Gaussian-distributed  random  vector.  The  above  model  is  well 
supported  by  experimental  data  provided  by  MIT  Lincoln  Laboratory  and  the  Jet  Propul¬ 
sion  Laboratory  in  the  form  of  polarimetric  measurements.  The  results  are  illustrated  by 
comparing  the  higher-order  normalized  intensity  moments  and  cumulative  density  func¬ 
tions  (CDF)  of  the  experimental  data  with  theoretical  results  of  the  K-distribution. 

Among  the  various  theoretical  models  applied  to  study  the  electromagnetic  wave 
scatterings  from  geophysical  terrain,  such  as  snow  and  ice,  the  radiative  transfer  theory 
has  drawn  intensive  attention  in  the  microwave  remote  sensing  society  during  the  past 
years.  In  most  of  the  scattering  models,  the  volume  scattering  and  the  surface  scattering 
effects  have  been  investigated  separately.  Recently,  there  has  been  a  growing  interest  in 
the  construction  of  composite  models  which  can  take  into  account  both  types  of  scattering. 
We  derived  the  first  order  iterative  solution  to  the  vector  radiative  transfer  equations  for  a 
two-layer  medium  with  a  diffuse  top  boundary  and  an  irregular  bottom  boundary  of  Gaus¬ 
sian  roughness.  The  Kirchhoff  approximation  and  the  geometrical  optics  approach  with 
shadowing  correction  are  used  in  formulating  the  boundary  conditions.  To  demonstrate 
the  utilities  of  the  theory,  randomly  oriented  spheroidal  discrete  scatterer  model  is  used  to 
calculate  the  backscattering  coefficients  from  soybean  field  in  different  growing  stages  and 
compared  to  the  experimental  measurements.  Good  agreement  has  been  achieved  for  both 
the  co-polarized  and  the  cross-polarized  data.  It  is  observed  that  the  presence  of  the  rough 
surface  can  significantly  enhance  the  backscattering  at  small  incident  angles  and  increase 
the  cross-polarized  returns.  The  polarization  signatures  calculated  based  on  the  Mueller 
matrix  show  a  straight  distortion  track  and  an  observable  pedestal.  Numerical  comparison 
to  the  backscattering  coefficients  calculated  by  using  planar  bottom  boundary  conditions 
with  or  without  the  incoherent  addition  of  the  rough  surface  effects  are  also  made. 
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The  concept  of  polarimetry  in  active  remote  sensing  is  extended  to  paissive  remote 
sensing.  The  potential  use  of  the  third  and  fourth  Stokes  parameters  U  and  V,  which 
play  an  important  role  in  polarimetric  active  remote  sensing,  is  demonstrated  for  passive 
remote  sensing.  It  is  shown  that,  by  the  use  of  the  reciprocity  principle,  the  polarimetric 
pairameters  of  passive  remote  sensing  can  be  obtained  through  the  solution  of  the  associated 
direct  scattering  problem.  In  particular,  the  full  polarimetric  information,  including  the 
corresponding  brightness  temperatures  of  U  and  V,  can  be  obtained  from  the  solution  of 
the  direct  scattering  problem  for  four  different  polarizations  of  the  incident  wave.  These 
ideas  are  applied  to  study  polarimetric  passive  remote  sensing  of  periodic  surfaces.  The 
solution  of  the  direct  scattering  problem  is  obtained  by  an  integral  equation  formulation. 
Incidence  on  a  penetrable,  lossy,  medium  is  considered.  Since  the  kernels  of  the  integral 
equations  are  the  periodic  Green’s  functions  and  their  normal  derivatives  on  the  surface, 
rapid  evaluation  of  the  slowly  convergent  series  associated  with  these  functions  is  observed 
to  be  critical  for  the  feasibility  of  the  method.  The  study  has  shown  that  the  brightness 
temperature  of  the  Stokes  parameter  U  can  be  significant  in  passive  remote  sensing.  Values 
as  high  as  50  K  are  observed  for  certain  configurations. 

To  demonstrate  the  use  of  polarimetry  in  peissive  remote  sensing  of  azimuthally 
asymmetric  features  on  a  terrain  surface,  an  experiment  was  designed  and  implemented.  A 
triangular  corrugation  pattern  was  made  on  the  sandy  soil  surface.  Polarimetric  brightness 
temperatures  are  measured  with  horizontal,  vertical,  and  45°  polarization  orientations  for 
various  observation  angles.  From  the  measured  temperatures,  absolute  values  as  high  as 
30-40  K  of  the  third  Stokes  brightness  temperatures  are  observed.  A  theoretical  analysis 
of  the  data  indicates  that  the  high  values  of  U  are  caused  by  the  azimuthal  asymmetry 
on  the  remotely  sensed  soil  surface.  It  is  also  observed  from  the  experiment  that  the 
brightness  temperatures  for  all  three  Stokes  parameters  vary  as  the  observation  direction 
varies  from  being  parallel  to  the  surface  row  structure  to  being  perpendicular  to  the  row 
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structure.  The  significant  implication  of  this  experiment  is  that  the  surface  asymmetry 
can  be  detected  with  a  measurement  oiU  at  a  single  azimuthal  angle. 
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Abstract 

To  interpret  the  polairimetric  active  remote  sensing  data  of  saline  ice,  v/e  have  developed  a 
random  medium  model  using  the  radiative  transfer  theory.  The  ice  layer  is  described  as  a 
host  ice  medium  embedded  with  randomly  distributed  inhomogeneities,  and  the  underly¬ 
ing  sea  water  is  considered  to  be  a  homogeneous  halfspace.  The  random  medium  model  is 
characterized  by  a  correlation  function  described  by  correlation  lengths  in  both  the  hori¬ 
zontal  and  vertical  directions.  Mtiltiple  scattering  effects  are  accounted  for  by  solving  the 
radiative  transfer  equations  numerically.  The  effects  of  random  roughness  at  the  air  -  ice, 
and  ice  -  water  interfaces  are  accoimted  for  by  modifying  the  boundary  conditions  in  the 
radiative  transfer  equations.  Analysis  of  the  model  for  reconstruction  of  sea  ice  parameters 
is  made.  An  optimization  approach  is  used  for  inversion.  The  discrepancy  between  the 
data  and  the  results  of  the  forward  model  is  minimized  by  changing  the  inversion  parame¬ 
ters  according  to  a  nonlinear  programming  scheme.  Reconstruction  of  correlation  lengths 
in  the  horizontal  and  vertical  dimensions  has  been  accomplished  using  the  polarimetric 
backscattering  coefficients  at  different  angles  of  incidence  as  input  data.  Effects  of  data 
diversity  and  noise  on  the  reconstruction  of  the  physical  parameters  of  sea  ice  from  the 
backscattering  coefficients  are  being  investigated. 


Polarimetric  Thermal  Emission  from  Randomly  Rough 

Ocean  Surfaces:  A  Numerical  Study 

J.  T.  Johnson,  R.  T.  Shin,  and  J.  A.  Kong 
Department  of  Electrical  Engineering  and  Computer  Science 
and  Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology,  Cambridge,  MA 

S.  H.  Yueh,  S.  V.  Nghiem,  and  R.  Kwok 
Jet  Propulsion  Laboratory 
California  Institute  of  Technology,  Pasadena,  CA 

Recent  theoretical  works  have  suggested  the  potential  of  passive  polarimetry  in  the 
remote  sensing  of  geophysical  media.  It  has  been  shown  that  the  third  Stokes  parameter 
U  of  the  thermal  emission  may  become  large  for  azimuthally  asymmetric  fields  of  obser¬ 
vation.  In  one  experiment,  values  of  U  of  as  hig^  as  40  K  were  measured  from  a  periodic 
triangular  soil  surface  at  10  GHz.  In  order  to  investigate  the  potential  applicability  of 
passive  polarimetry  to  the  remote  sensing  of  ocean  surface,  a  numerical  study  of  the  po¬ 
larimetric  thermal  emission  from  randomly  rough  ocean  surfaces  was  performed.  A  Monte 
Carlo  technique  utilizing  an  exact  method  for  calculating  thermal  emission  was  chosen  for 
the  study  to  avoid  any  of  the  limitations  of  the  commonly  used  approximate  methods  in 
rough  surface  scattering. 

In  this  Monte  Carlo  technique,  a  set  of  finite  rough  siirface  profiles  in  two  dimensions 
with  desired  statistics  was  generated.  Each  fimte  surface  was  then  extended  periodically 
to  create  a  set  of  infinite  rough  surfaces.  The  polarimetric  thermal  enussion  from  each 
surface  of  the  set  was  then  calculated  tiring  both  the  extended  botmdary  condition  method 
and  the  method  of  moments.  Finally,  the  results  from  the  set  were  then  averaged  to  obtain 
the  estimate  of  the  polarimetric  brightness  temperatures  for  the  given  surface  statistics. 
The  surface  statistics  chosen  were  intended  to  model  a  wind  perturbed  ocean  surface  in 
the  X  to  band  microwave  re^on.  For  this  purpose,  a  power  law  surface  spectrum 
corresponding  to  the  capillary  wave  portion  of  the  ocean  spectrum  was  chosen  with  rms 
surface  heights  corresponding  to  wind  speeds  of  up  to  10  m/s.  The  effects  of  var^dng  the 
azimuthal  and  polar  looking  angles,  rms  surface  bright,  ocean  permittivity,  exponent  of 
the  power  law  spectrum,  and  the  necessary  high  and  low  frequency  cutoffs  of  the  spectrum 
were  investigated  in  this  numerical  study. 

The  results  of  the  study  indicate  that  the  U  parameter  is  sensitive  to  the  aizimuthal 
angle  between  the  surface  periodicity  and  the  looking  angle  and  to  the  rms  height  of 
the  surface,  and  that  the  U  parameter  is  fairly  insensitive  to  variations  in  polar  angle, 
permittivity,  surface  power  law  spectrum,  and  surface  spectrum  high  frequency  cutoff. 
These  properties  pre  further  strength  to  the  idea  of  uring  the  U  parameter  to  detect  wind 
direction  over  the  ocean. 
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Abstract  -  This  paper  presents  an  anisotropic  layer  model  with  ellipsoidal  scatterers  for 
applications  to  polarimetric  remote  sensing  of  geophysical  media  at  microwave  frequencies. 
The  scattering  configuration  includes  an  isotropic  layer  covering  an  anisotropic  layer  above 
a  homogeneous  half  space.  Tb<*  isotropic  layer  consists  of  randomly  oriented  spheroids. 
The  anisotropic  layer  contain^  ellipsoidal  scatterers  with  preferential  vertical  alignment 
and  random  azimuthal  orientations.  Effective  permittivities  of  the  scattering  media  are 
calculated  with  the  strong  permittivity  theory  extended  to  account  for  the  non-spherical 
shapes  and  the  scatterer  orientation  distributions.  Based  on  the  analytic  wave  theory, 
the  dyadic  Green’s  function  for  layer  media  is  used  to  derive  polarimetric  backscattering 
coefficients  under  the  distorted  Bom  approximation.  The  ellipsoidal  shape  of  the  scatterers 
gives  rise  to  non-zero  cross-polarized  returns  from  the  untilted  anisotropic  medium  even 
in  the  first  order  approximation.  The  effect  of  rough  interface  is  estimated  by  incoherent 
addition  method.  Theoretical  results  and  experimental  data  are  matched  at  9  GHz  for 
thick  first-year  sea  ice  with  bare  surface  and  with  snow  cover  at  Point  Barrow,  Alaska. 
The  model  is  then  used  to  study  the  sensitivity  of  the  scattering  coefficients  with  respect 
to  correlation  lengths  representing  the  geometry  of  brine  inclusions  in  the  model  for  sea 
ice.  Finally,  polarimetric  signatvires  of  bare  and  snow-covered  sea  ice  are  also  simulated 
based  on  the  model  to  investigate  effects  of  different  scattering  mechanisms. 
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1.  INTRODUCTION 

Technology  in  remote  sensing  has  been  advanced  considerably  especially  for  airborne 
and  space  borne  radar  with  multi-frequency  and  multi-polarization  capabilities.  Theoreti¬ 
cal  models  have  been  developed  to  interpret  multi-frequency  polarimetric  data  for  remote 
sensing  of  geophysical  media.  Several  approaches  including  radiative  transfer,  modified  ra¬ 
diative  transfer,  and  ainalytic  wave  theory  have  been  considered  in  model  developments  [1]. 
The  analytical  wave  theory,  while  remains  mathematic2Jly  tractable  for  some  complexity  of 
the  media,  preserves  the  phase  information  which  is  appropriate  for  the  calcidation  of  po¬ 
larimetric  scattering  coefficients  [2].  Under  the  Born  approximation  for  sparse  and  tenuous 
media,  conventional  scattering  coefficients  have  been  derived  for  isotropic  multi-layer  [3], 
anisotropic  layer  [4],  and  isotropic-anisotropic  layer  [5]  configurations.  Fully  polarimetric 
scattering  coefficients  have  also  been  calculated  for  isotropic  [6]  and  anisotropic  [7]  layer 
configurations.  In  these  models,  the  scatterers  are  spherical  or  spheroidal  and  the  cross- 
polaxized  returns  come  from  second  (or  higher)  order  [8]  or  the  tilted  anisotropy  [4].  For 
denser  or  less  tenuous  media,  the  distorted  Born  approximation  hzts  been  applied  [9-11]. 
This  approximation  considers  dissipation  loss,  scattering  loss,  and  the  modification  of  wave 
speed  due  to  the  scatterers;  thus,  multiple  scattering  has  been  included  to  some  extent. 
For  media  with  strong  permittivity  fluctuations,  the  strong  permittivity  fluctuation  theory 
is  used  in  conjunction  with  the  distorted  Born  approximation  [12].  In  this  case,  conven¬ 
tional  backscattering  coefficients  have  been  computed  for  an  isotropic  half-space  [13],  an 
anisotropic  half-space  [14],  and  a  configuration  for  an  isotropic  layer  of  spherical  scatterers 
above  an  anisotropic  layer  of  aligned  spheroids  [15]. 

In  this  paper,  the  isotropic-anisotropic  layer  configuration  containing  non-spherical 
scatterers  is  considered  to  obtain  fxiUy  polarimetric  backscattering  coefficients  for  geophys¬ 
ical  media.  The  scattering  configuration  is  illustrated  in  Figure  1  where  6oi  is  the  incident 
angle.  The  covering  isotropic  layer  is  composed  of  randomly  oriented  spheroids  as  described 
in  reference  [16].  In  the  anisotropic  layer,  the  scatterers  are  modeled  with  sm  ellipsoidal 
correlation  function  with  the  orientation  characterized  by  a  probability  density  function 
of  the  Eulerian  rotation  angles.  The  orientation  of  the  ellipsoids  is  vertically  aligned  and 
aaimuthally  random.  The  strong  permittivity  fluctuation  theory  is  extended  to  calculate 
the  effective  permittivities  and  the  distorted  Born  approximation  is  applied  to  derive  the 
polarimetric  scattering  coefficients.  Theoretical  results  are  compsired  with  measured  data 
for  bare  and  snow-covered  sea  ice.  The  effects  of  rough  interfaces  aire  also  estimated  by 
the  incoherent  addition  approach  [4].  Variations  of  polarimetric  backscatter  are  studied 
for  various  correlation  lengths  and  polarization  signatures  of  sea  ice  are  simulated  and 
discussed  for  different  scattering  mechanisms.  Following  this  plan,  the  paper  consists  of 
5  sections.  Section  2  is  for  effective  permittivities,  Section  3  for  polarimetric  scattering 
coefficients.  Section  4  for  results  and  discussion,  and  finally  Section  5  for  the  summary. 
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2.  EFFECTIVE  PERMITTIVITY 
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In  the  isotropic  scattering  medium  of  Region  1  in  Figure  1,  the  embedded  scatterers 
are  modeled  as  spheroids  and  the  effective  permittivity  has  been  derived  tmder  the  strong 
fluctuation  theory  and  reported  in  [16,17].  In  the  anisotropic  medium,  the  ellipsoidal 
scatterers  have  preferentied  alignment  resulting  in  the  effective  zinisotropy.  The  effective 
permittivity  tensor  of  the  scattering  anisotropic  medium  will  be  calculated  in  this  section 
for  the  ellipsoids  with  vertical  alignment  and  random  azimuthal  orientations. 

In  an  inhomogeneous  medium  such  as  sea  ice,  sea  water  is  trapped  in  an  ice  medium 
in  form  of  brine  inclusions  which  are  usually  ellipsoidal.  The  ice  tends  to  grow  verti- 
caJly  downward  rendering  the  ellipsoidal  inclusions  aligned  preferentially  in  the  vertical 
direction.  In  the  absence  of  sea  currents,  the  crystallographic  c-axes  are  random  in  the 
horizontal  plane  as  seen  in  Figure  2,  which  depicts  a  horizontal  thin  section  of  sea  ice. 
Consequently,  the  minor  axes  of  the  ellipsoids  have  random  orientations  parallel  to  the 
horizontal  plan.  In  this  model,  a  correlation  function  corresponding  locally  to  a  scatterer 
is  used  in  the  derivation  of  the  effective  permittivity  with  the  strong  permittivity  fluctu¬ 
ation  theory,  extended  to  account  for  the  orientations  of  the  ellipsoidal  scatterers.  When 
the  average  process  is  performed  over  orientation  angle  <f>f  shown  in  Figure  3,  the  effec¬ 
tive  permittivity  is  an  untilted  uniaxial  tensor  with  vertical  optic  axis,  which  effectively 
manifests  the  azimuthal  symmetry  of  the  inhomogeneous  medium. 

Let  €1,  be  the  permittivity  of  the  host  medixim  2uid  e,  be  the  permittivity  of  the  em¬ 
bedded  ellipsoidal  scatterers  occupying  a  total  fractional  volume  of  The  subscript  for 
the  anisotropic  medium  (Region  2)  is  omitted  in  this  section  for  convenience  since  the  fol¬ 
lowing  derivation  is  for  Region  2  only.  Similar  to  the  method  in  [12],  auxiliary  permittivity 
Ig  =  diag[€gp,€gp,€gz]  Is  iutroduced  into  the  wave  equation  for  latter  consideration  of  the 
singularity  in  anisotropic  dyadic  Green’s  function  The  singularity  is  accounted 

for  by  decomposing  Gg{f,r*J^  =  Gg(f  —  r*)  into  a  principal  value  part  and  a  Dirac  delta 
part  with  dyadic  coefficient  S  which,  in  this  case,  is  a  diagonal  tensor  with  three  distinctive 
diagonal  elements  Syi,  and  Sz>  in  the  local  coordinates  (x',y',z').  These  coordinates 
are  related  to  the  global  coordinates  {x,y,z)  by  the  Eulerian  rotation  tensor 

cos<l>f  sm<f>f  O' 

T  =  —sin4>f  cos<f>f  0  (1) 

0  0  1. 

The  effective  permittivity  of  the  inhomogeneous  medium  is  composed  of  a  quasi-static 
part  and  a  scattering-effect  part  which  accounts  for  the  attenuation  and  the  modiflcation 
in  the  wave  speed  due  to  the  inhomogeneities 

«e//  =  +  €o  -  ?e//  ‘  •  leff 


(2) 
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where  auxiliary  permittivity  €g  and  dyadic  coefficient  which  is  the  ensemble  average 
of  5,  Me  determined  by  the  condition  of  secular-term  elimination.  The  effective  dyadic 
scatterer  under  the  low-frequency  approximation,  is  given  by 


_  *2ir  *»Vi*  f  foo  _  _  N  I 

=  l  [G,(E')],, $«(?)+ [5]„l  (3) 

in  which  p(^/)  is  the  probability  density  function  of  orientation  angle  is  the 

variance  defined  in  [2],  fco  is  the  free-space  wave  number,  and  is  the  Fourier  transform  of 
the  normalized  local  correlation  function.  The  anisotropic  Green’s  function  Gg  [12],  which 
is  invariant  under  the  azimuthal  Eulerian  rotation  (1),  is  expressed  in  the  k'  domain  as 


G,{k  )  = 


(k'l  +  k'l)D.(k ) 


r  k'l 
-KK 
L  0 


{k'l  +  k'l)D.{k  ) 


r 


KK 


KK 

I.*2 


-KK 


0 

0 

^Mk’l  +  k'l) 


+ 


k  k 


kIMk  ) 


(4a) 


where  the  quantities  Do(k  ),  De(k  ),  and  k^^  are  in  the  following  equations  with 
angular  frequency  u  and  permeability  fio 


D,(t)  =  k'l+k'l  +  k'l-kl, 

(4b) 

D.(^)  =  k'l  +  '^^(k']:+k-l-kl) 

(4c) 

^BP  -  and  fcj,  =  (v^noCgt 

(4d) 

As  in  (2),  the  effective  permittivity  has  been  approximated  by  truncating  the  series  in  the 
renormalization  method.  The  validity  condition  for  the  approximation  is 

The  ellipsoid^J  scatterer  is  described  with  a  normalized  local  correlation  function  of 
the  form  _ 

Bj(r')=«p(-^J  +  ^  +  ^)  (6a) 

with  correlation  length  £*» ,  £yi ,  and  in  the  local  coordinates  corresponding  to  the  minor, 
the  meridian,  and  the  major  axes  of  the  scatterer.  In  this  model,  the  correlation  lengths 
are  related  to  the  effective  size  and  shape  of  the  scatterers.  This  local  correlation  function 
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can  be  reduced  to  spheroidal  shape  for  two  equal  correlation  lengths  and  to  spherical  shape 
for  three  identical  correlation  lengths.  Applying  the  Fourier  transform  on  (5a)  yields 

For  random  horizontal  orientations  with  no  preference  in  azimuthal  direction,  the  proba¬ 
bility  density  function  of  orientation  is  simply 


P(^/)  =  l/(27r)  (6) 

To  calculate  the  effective  permittivity  according  to  (2),  and  (5)  need  be  deter¬ 
mined.  Due  to  the  global  azimuthal  syirimetry,  auxiliary  permittivity  f ^  in  the  coordinate 
(x,y,  z)  is  uniaxial  as  afore  indicated  and  the  elements  in  Ig  are  subjected  to  the  condition 
(^)  =  0  such  that 


f 


0 

0 


0 

U'  +  iv' 
0 


)  =  0 


(7) 


Local  quantities  and  ^gi  in  (7)  are  related  to  the  elements  of  dyadic  coefficient 

5  =  dia^[5a,», 5j,», 5z*]  by 


^  f  \  ^  ^SP 

eo+SAe- egg) 

(8a) 

e,+Sy>{e-egp) 

(8b) 

(8c) 

where  c  can  takes  on  the  value  of  e,  in  a  scatterer  or  in  the  background  medium.  From 
(7)  and  (8),  Cgp  and  Cj*  can  be  written  as,  respectively. 


_  ,  ^gp  2eo  d-  (Sg*  +  5'y>)(e,  Cgp) 

I- f,  2co  +  iSg.  +Sg>)iei-  egp) 

gP  +  5«>(gt  —  £gp)  ep  +  Sy>{€l,  —  Cgp) 

Cp  +  Sx'itg  —  ^gp)  Cp  +  Sy>{^€g  — 
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—  €6  + 


1-/- 


gp  d-  Sz>{€i  —  Cgg) 

€p  5'*»(c»  ^gz) 


(9a) 

(9b) 
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The  average  dyadic  coefficient  in  the  global  coordinates  is  obtained  by  the  averaging 
integration  over  the  probability  density  function  of  orientation 

_  /•2ir  ^  ^  = 

(S)  =  /  d(f> f  p{<f> f)  T  •  0  Sgi  0  'T 

do  n  n  c 


0  0 


'  Sx>  +  Sy> 

0 

0  ■ 

\Sp 

0 

0  ■ 

0 

5**  +  Sy> 

0 

=  0 

0 

0 

0 

2Sz'. 

0 

0 

5.. 

The  coefficients  Sx>,  Sy‘,  and  Sz'  are  derived  from  the  seculair  elimination  condition  [12] 
which  renders 


Jo  2iregpay/a 

s,.  =  [(1  +o)lan->  VS] 

Ja  2ir«j,aV“ 

5s'  =  (  d4>  _  tan"^ 


(lla) 

(llb) 

(llc) 


where  the  values  of  both  the  squzure  root  and  the  inverse  tangent  are  chosen  on  the  principal 
Riemann  sheets  with  branch  cuts  on  the  negative  real  axis.  In  (11),  the  integrations  over 
<}>  can  be  carried  out  numerically  and  quantities  a,  7^,  and  7y  axe  defined  as 


2  1 

a  =  aqf  —  1 ,  a  = 

^gp 


_  1  I  cos^  ^  sin^  <f>  \ 
\  ( cos^  <f>  sin^  <f>\ 


_  1  /  cos^ 


<f)  sin^ ' 
^y' 


(12a) 


(12b) 


(12c) 


(12d) 


Also  due  to  the  azimuthal  symmetry  of  the  scattering  medium,  the  effective  scatterer 
tensor  ^eff  has  the  uniaxial  form 

_  6//P  0  0 

(eff  =  0  Ceffp  0  (13) 

L  0  0  ^e//z  . 
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which  is  obtained  by  substituting  the  Green’s  function,  the  correlation  function,  and  the 
probability  density  function  of  orientation  into  (3).  The  integrations  are  caxried  out  with 
the  procedure  in  [2,17].  Followed  is  how  to  compute  : 


^e//p  =  2  +  ^y')] 

^e//z  =  (Iz<  +  St') 


where  variance  ,  8^y> ,  and  S^t'  are  respectively  determined  by 


— 


S^y>  = 


8fz>  = 


€6  ^gp 

(1-/.)+ 

C,  Cgp 

€0+<S'*»(c6  — €jp) 

£o+‘S'*'(Cj  — Cjp). 

Cfc  ^gp 

(1-/.)+ 

^t  —  ^gp 

.^0+Sy>(£b-egp), 

.Co +5'-'(e,  —  Cpp) 

€b  Cgz 

(1-/.)+ 

~  ^gz 

^€Q+Sz'{eb  —  £gz). 

£o-\-Sz'{it—^gz)  _ 

fz 


f. 


f. 


(14a) 

(14b) 


(15a) 

(15b) 

(15c) 


As  observed  from  (14),  expressions  for  /*<,  Iy>^  and  It'  are  necessary  to  complete  the 
derivation  of  the  anisotropic  effective  permittivity.  The  result  for  Iz>  is 


A.=  rd4,^(xi+xi) 


(16») 


(l_  tail-* 


*  2a2  [  T?e  V2 


J5  = 


t?o(a+2)-(6+ai/J,,) 

2a? 

t?o 

(^--tan 


vrj\ 


=  C=«7\^. 

+  C 


b  = 


,  i?„  =  6-  l,  i?^  =  6  +  C 


For  /*' ,  the  result  is  expressed  as  follows 

k]i 


Ix>  =  /  d<f>—  sin*^  +  ilij*  cos'^<f>  I^) 

Jo  ^ 


/;=«  (ii-z5_2j)  +  l(i,+zi_l) 


(16b) 

(16c) 

(16d) 

(16e) 

(17a) 

(17b) 

(17c) 
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25'  =  2Tla  =  l.  2?  =  2Jlc<  =  l,  2?=2|la  =  l  UW 


25'  =  2Tla  =  l.  2?  =  2Jlc<  =  l,  2?=2|la  =  l  UW 

From  the  symmetry,  Explicitly,  the  result  for  /yi  is 

-  f  d4>—  {ll.'rt  cos*^  1%  +  sin^^  1%)  (18) 

Substituting  fj,  and  in  (2)  yields  tbe  uniaxial  effective  permittivity  tensor 
Ce//  whose  lateral  and  vertical  elements  are,  respectively 

Ce//P  =  (■gp  +  ^oieffp/i)-  ~  Sp^effp)  (19a) 

^effz  —  ^gz  +  €0^e//*/(l  ~  ^z^effz)  (191>) 


As  seen  from  the  above  expressions,  effective  permittivity  I^ff  is  anisotropic  with  optic 
axis  in  the  vertical  direction.  In  the  next  section,  the  anisotropic  effective  permittivity  is 
used  in  the  derivation  of  the  polaximetric  backscattering  coefRcients  under  the  distorted 
Born  approximation  with  dyadic  Green’s  function  (DGF)  for  the  layer  configuration. 

3.  SCATTERING  COEFFICIENTS 

Consider  the  scattering  configuration  in  Figure  1.  Region  0  is  the  upper  half-space 
with  permittivity  cq  •  Region  1  of  thickness  di  is  an  isotropic  medium  consisted  of  randomly 
oriented  spheroidal  scatterers  of  permittivity  e«i  in  a  background  medium  of  permittivity 
£{,1 .  Region  2  of  thickness  dj  is  an  anisotropic  medium  composed  of  ellipsoidal  scatterers  of 
permittivity  c,2  embedded  in  a  host  medium  of  permittivity  ei2-  The  ellipsoids  axe  oriented 
preferentially  in  the  vertical  direction  and  randomly  in  azimuthal  directions  as  described 
in  the  last  section.  Region  3  is  the  underlying  homogeneous  half-space  of  permittivity 
C3.  The  effective  pernoittivities  of  the  scattering  regions  are  calculated  as  shown  in  the 
Izist  section.  The  distorted  Born  approximation  is  applied  to  derive  the  complete  set  of 
polarimetric  backscattering  coefficients  constituting  the  covariance  matrix  and  the  Mueller 
matrix  characterizing  the  polarimetric  scattering  properties  of  the  layer  media. 

Polairimetric  backscattering  coefficients  have  been  defined  with  ensemble  averages  of 
scattered  fields  [2].  The  averages  are  calculated  with  spatial  integrations  over  products  of 
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the  DGFs,  the  mean  fields,  and  the  correlation  functions  as  follows 

(Eoa(r)  •  ;^,(r))  =  f  dil;f  f 

i  do  Jo 

JVi  JVi 


+  ^0  /  d<j>f  P2{4>f)  dT2  I  dr^C(2jklm(r2,f^;<f>f) 

ijXl.m  -/o  Jy^ 

•  [(<?02tj(^,»^2))(-F2l!(F2))]  •  [(Go2t/(r,f|))(F2m(r|))] 


(20) 


where  V\  is  the  volume  occupied  by  Region  1  and  V2  by  Region  2.  The  DGFs  and  the 
mean  fields  have  been  obtained  in  [2];  the  correlation  functions  need  to  be  specified  next. 
The  integrations  are  then  effectuated  to  derive  the  scattering  coefficients. 

For  Region  1,  the  correlation  functions  C7*s  in  (20)  zire  defined  in  [16]  for  the  spheroids. 
In  Region  2,  the  correlation  functions  in  the  global  spatial  domain  are 

G(2iifelm(r2,rf,^/)  =  (6ifc(»^2)^2<m(»t)|<^/(»^2))  (21) 

To  facilitate  the  integration  of  (20),  Fourier  transforms  of  the  correlation  fimctions  are 
introduced  for  the  statistically  homogeneous  scattering  media  under  consideration 

=  r  -  *1)  (22) 

j  —00 

which  is  expressed  in  the  global  coordinate  system  (x,y,z)  and  related  by  the  Eulerian 
rotation  transformation  T  (1)  to  the  following  non-zero  correlations  in  the  local  coordinate 
system  (x',y',z') 


^2z'*'(^  )  =  $2*'«»*'*'(/^  )  =  ) 

(23a) 

^2x'y'ifi  )  =  ^2*'«*y'y'(^  )  =  hx'y'^tZ^P  ) 

(23b) 

^2x'z'{^  )  =  ^2x’x’x’z'{P  )  =  ^2*'*'^(2(j9  ) 

(23c) 

$2y'i'(/3  )  =  ^2y>y>x'x>{P  )  =  ^2y'i'^{2(/3  ) 

(23d) 

$2y'y‘(^  )  =  ^2y»y'y'y»(/3  )  =  hy’y>^i2{^  ) 

(23e) 

^2y'*'(^  )  =  ^2y'y» *'«'(/?  )  =  hy>z’^i2{fi  ) 

(23f) 

^2z'x'{P  )  =  ^2z'z'x»x'{P  )  =  ^2z'e»${2(^  ) 

(23g) 

$2*'y'(^  )  =  ^2*'*'y»y'(^  )  =  ^2*'y'^(2(/3  ) 

(23h) 

^2*'z'(^  )  =  ^2z'z'z'z'(^  )  =  ^2z’z'^(,2{P  ) 

(23i) 
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where  )  is  given  in  (5b),  which  is  the  Fourier  transform  of  (5a),  and  the  variance 

^2*8  in  (23)  are 


S2x‘z'  = 


82y>y>  — 


82z>z‘  = 


^b2—^g2p 

Co  ■\-S2x>{^h2—^g2p) 

€62  ~  ^g2p 

Co  +  S2y>  (C62  —  ^g2p) 

€62  — €g2z 

Co  +  5'22»  (C62  —  ^g2z) 


(1  — /*2)  + 

€,2  —  (■g2p 

Co  +  i52*'(c«2  —  Cs2p) 

(1— /•2)  + 

£.2—£g2p 

Co  +  S2y>  (c,2  —  ^g2p) 

\l-f.2)  + 

C*2  —  Cp2z 

Co  +  S2z'  (^.2  —  ^g2z) 

f.2 


f»2 


f.2 


82x'y'  = 


+ 


^2x'z’  = 


+ 

82y>z'  = 

+ 


^b2—^g2p 

^62  ^p2p 

Co  +  52*'  (€62  —  ^g2p) . 

Co  +52y'  (C62  —  ^g2p) . 

€j2  — Cj2p 

Co  +  52*' (c,2  ~^g2p). 

Co  +  S2y>  (c,2  —  ^g2p). 

^b2—^g2p 

^62— «p2z 

Co  +  52*'  (C62  —  ^g2p) . 

Co  +52*'(c62  —  Cp2z) . 

€.2  —  Cy2p 

C/2  —  ^g2z 

Co  +  52*'  (€,2  —  Cj2p) . 

_€o  +  52«'(c,2  — €y2*). 

€b2—€g2p 

CM— ep2z 

Co  +  52  v'  (C62  —  Cp2p) . 

Co  +  52z'  (cm  —  Cy2z) . 

€*2  — €ff2p 

C»2  Cy2z 

Co  +  S2y>  (€,2  —  €p2p) . 

Co  +52z'(c,2  —  Cp2z) . 

(1-/-2) 


ft2  =  ^2y'x' 


(1-/.2) 

♦ 

f.2  =  ^2z'x' 


(1-/.2) 


f.2  =  82glyl 


(24a) 

(24b) 

(24c) 


(24d) 


(24e) 


(24f) 


Based  on  the  invariant  property  of  the  Fourier  transform  under  the  rotation  trsmsformation, 
spectral  density  ^2jklm{P)  in  the  global  coordinates  can  functionally  be  related  to  those 
given  in  (24)  with 


$(2(/?)  =  ^i2{Px  =  P*  COS^/  +  Py 

P'y  =  -Px  sin  ^/  +  Py  cos  4>fi 

P'z=P.)  (25) 

In  the  global  coordinates,  the  rotation  transformation  together  with  the  above  invariant 
property  cast  the  anisotropic  spectral  densities  into  the  form 


^2jklm{P)  =  82jklm^i2{P) 


(26) 


V 


where  variance  62jklm  are  dependent  on  the  Euierian  angle  as 


^2jklm  ~  {,^2z'x‘Txxjk  d"  ^2y‘x'1‘yyjk  "!■  ^2x'x'Tzzjk'}Txxlm 
"H  {^2x'y’Txxjk  *H  ^2y‘y'Tyyjk  "h  ^2z'y'Tzzjk')  '^yylm 
+  {Six' z'Txxjk  "h  Siy! z'Tyyjk  +  Siz' z'^ZZ^k)  ^zzlm  (27) 

The  Euierian  rotation  has  been  applied  to  arrived  at  (27)  where  T’s  are  elements  of  the 
following  transforming  tensors 

cos^<f>f  sin^/cos^/  O' 

Txx  =  sin  cos  sm^<f>f  0  (28a) 

0  0  oj 

six?  4>f  — sin^/cos^/  O' 

Tyy  =  —  sin^/cos^/  CO?  <f>f  0  (28b) 

0  0  0. 


(28a) 


(28b) 


_  _  _  _  0  0  0 
Tzz—I  'I'xx  ^yy  ~  0  0  0 

0  0  1 


(28c) 


By  substituting  the  above  correlation  functions,  the  dyadic  Green’s  functions  of  the 
layer  medium,  and  the  mean  fields  into  (20),  the  correlations  of  the  scattered  field  can  now 
be  found.  Rearran^ng  all  the  coefficients  from  the  DGF  and  mean  fields,  the  scattering 
coefficients  can  be  written  as 


-1,1  x,y,z 


=  irkt  j  dij^f  j  d<f>f  pi{iPf,<f>f)  ^  ^ 

®  **  a,bfC,dj,k,l,m 


ab  n-ahed 

lfir,jk^lvK,lm-^ljklm 


ou,od 

,2-K  eu,ed  x,y,z 

The  scattering  contribution  of  the  first  term  in  (29)  &om  Region  1  due  to  the  spheroids 
has  been  obtained  in  [16]  and  the  solution  consists  of  16  terms  for  downgoing  and  upgoing 
incident  and  scattered  waves.  Region  2  contains  scatterers  with  preferentizd  vertical  align¬ 
ment,  which  effectively  renders  the  medium  anisotropic.  The  ordinary  and  extraordinary 
waves  going  down  and  up  give  rise  to  216  terms  in  the  second  term  of  (29).  All  coefficient 
'J'a’s  have  been  derived  and  given  in  [2].  For  the  result  is 

_  ,262jklmS2x'i2y'^2z'  T  e-*(«F«  “  •^r,)d2 

-  *  .A 


'’2jklm 


irl\z. 


e” *(«pq  —  f 
{lipq  ^2'?{^pq 


gi(/Cr«  — 

(^rj  “  ^2?{^r»  ~  1^2  )^(^»“*  — 


*)*(''?«  “  ^rs) 


—  P2(/Ca)  —  fi2('C2) 


(30a) 
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where  'P2(«2)  and  Q2('C2  )>  resulted  from  complex  integrations,  are  determined  by,  respec¬ 
tively. 


^2(«2)  = 


+ 


t(d2  -  di)c“*('‘2  -  /Cr.)dlgt(/C2  “  ACp,)d2 
(2tIm/C2)*(«2  —  '«pj)(«2  —  «rj) 

■g— t(/Cpj  —  /Cr«)di  g“*(^P9  ^r*')^2 

(2tlm«2)^(«2  —  'Cp9)('C2  —  Krs) 

e-*(«2  -  «r#)dlgt(/C2  -  Kp,)d2  ' 
(2tIm/C2)*('C2  —  «P9)(^2  —  'Cm). 


1 

tim  /C2 


+ 


«2  —  >ipq  'C2  —  'Cr» 


Q2('C2)  — 


_  i(d2  -  di)c*(''2  -  '«P9)‘^lc-*'('«2  -  >^T,)d2 


+ 


(2tImK*)2(/c*-/Cp,)(/c*-/c„) 

C*(^2  ~ “ 'Cr»)d2 
(2llm /C2  )^('C2  “  '®i>9)('^2  ~  'Crj) 


ilm  /ci 


♦  + 


'^2  'Cp9 


+ 


/Cj  ~  'C 


•r» 


In  (30),  the  pole  K2  is  computed  with  the  following  equations 


«2  =  i^2l'^'^  +  klAx'  +  ^lv^y' 

=  kpi  cos{4>i  -  4)f)  -  kp,  cos{<f>, 
k^y  =  kpi  sin(^{  ^/)  kpt  sin(^f  —  ^/) 


(30b) 


(30c) 


(31a) 

(31b) 

(31c) 


where  <l>i  and  <f>g  are  aizimuthal  angles  of  incident  and  scattered  waves,  respectively.  The 
remedning  integrations  over  Eulerian  angles  in  (29)  can  be  carried  out  numerically.  When 
low-frequency  condition  is  valid,  the  integrations  can  be  done  analytically.  It  should  be 
noted  that  (29)  expresses  the  scattering  coefficients  in  the  scattered  basis  which  can  be 
transformed  to  the  incident  basis  by  changing  the  sign  of  ahTVK  and  <r„rhK  [2].  The  scat¬ 
tering  effect  of  the  rough  boundaries  at  the  mediiim  interfaces  can  also  be  estimated  by 
incoherently  adding  the  total  contribution  from  the  rough  surface  scattering  with  consider¬ 
ation  of  the  propagation  loss  in  the  calctdations  of  scattering  coefficients.  Compared  with 
the  model  in  [2]  for  the  case  of  vertically  oriented  spheroids,  the  new  model  can  give  higher 
co-polarized  backscattering  coefficients  due  to  higher  total  cross  section  of  the  ellipsoids  for 
the  same  fractional  volume.  Another  difference  is  that  the  new  model  provides  non-zero 
cross-polarized  backscattering  in  the  first-order  distorted  Bom  term  whereas  the  former 
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can  only  account  for  the  cross-polarized  return  in  the  second  or  higher  order  term  when 
the  spheroids  are  vertically  oriented. 

4.  RESULTS  AND  DISCUSSION 
a.  Data  Comparisons 

In  this  subsection,  theoretical  results  are  compared  with  experimental  measurements 
for  backscattering  coefficients  at  9  GHz  as  a  function  of  incident  angle.  The  radar  backscat- 
ter  data  were  collected  for  the  thick  first-year  sea  ice  near  Point  Barrow,  Alaska  [86].  For 
thick  first-year  sea  ice  with  bare  surface,  the  scattering  configuration  is  illustrated  in  Figure 
4  representing  a  sea  ice  layer  over  sea  water.  The  ice  layer  was  1.65-m  thick  and  contained 
brine  inclusions.  From  ground-truth  data,  it  has  been  inferred  [19]  that  the  brine  had  a 
permittivity  of  =  (38.0  -f  41.0)£o  &Qd  occupied  a  fractional  volume  of  f,i  =  4.5%.  The 
background  ice  permittivity  was  estimated  to  be  =  (3.15-i-t0.002)£o  and  the  permittiv¬ 
ity  of  sea  water  £3  =  (45.0  +  i40.0)£o.  With  correlation  length  =  0.70  mm,  =  0.25 
mm,  and  iiz>  =  1*20  mm,  the  theoretical  results  match  well  with  the  experimental  data,  as 
shown  in  Figure  5,  for  co-polarized  backscattering  coefficient  akk  at  large  incident 

angles,  and  cross-polarized  backscattering  coefficient  over  the  range  of  incident  an^es. 

The  disagreement  in  the  co-polarized  backscattering  coefficients  at  small  incident 
angles  is  due  to  the  scattering  from  rough  surface.  In  the  configuration  of  Figure  5,  the 
interfaces  have  been  assumed  to  be  smooth  and  the  rough  surface  scattering  has  been 
ignored.  The  effect  of  the  roughness  at  the  lower  interface  of  the  thick  first-year  sea  ice 
is  negli^ble  since  the  wave  at  14  GHz  is  attenuated  before  the  sea  water  is  reached.  To 
account  for  the  rough  surface  effect  at  the  top  interface,  a  Gaussian  roughness  is  now 
considered  with  standard  deviation  <r,  =  0.6  Tnm  and  correlation  length  I,  =  1.5  cm  as 
depicted  in  Figtire  6,  where  all  other  physical  parameters  are  the  same  as  in  Figure  4. 
The  rough  surface  contribution,  calculated  with  the  small  perttixbation  method  (SPM)  [1], 
is  incoherently  added  to  the  volume  scattering.  The  comparison  between  theoretical  and 
experimental  results  for  the  co-polarized  returns  is  improved  at  the  low  incident  angles  as 
seen  in  Figure  7.  For  this  surface,  the  additional  contribution  to  the  cross-polarized  return 
is  small  and  is  actually  ignored  in  the  SPM  model  applied  in  this  case.  The  cross-polarized 
return  therefore  remains  unchanged. 

The  configuration  for  snow-covered  sea  ice  is  shown  in  Figure  8.  The  snow  is  10- 
cm  thick  according  to  the  average  of  the  thickness  range  reported  in  [18].  In  the  snow 
layer,  the  oblate  spheroidal  shape  is  assumed  for  the  ice  grains  with  correlation  lengths 
=  ^x/6  =  0.15  mm.  The  physical  parameters  for  the  sea  ice  layer  and  the  sea  water 
are  kept  unchanged.  Both  the  top  and  the  nuddle  interfaces  are  rough  with  stsmdard 
deviation  a,  —  0.1  cm  and  correlation  length  s:  1.5  cm.  The  theoretical  result  compare 
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well  with  the  measurements  from  the  sea  ice  layer  with  snow  cover  2is  seen  in  Figure  9  for 
the  conventional  backscattering  coefficients.  The  theoretical  cross-polaxized  returns  are, 
however,  lower  tham  the  measurements.  For  a  medium  with  high  fractional  volume  such 
as  snow,  a  significant  part  of  the  total  contribution  to  Chv  can  come  from  higher  order 
scattering  which  is  ignored  in  the  present  model.  Regarding  the  scattering  mechanisms, 
volume  scattering  is  dominant  at  large  incident  angles.  Rough  surface  effect  is  still  strong 
at  small  incident  angles  with  more  contribution  from  the  snow-ice  interface  due  to  the 
larger  permittivity  contrast.  The  oscillations  observed  in  the  curves  in  Figure  9  sire  caused 
by  the  boundary  effect  of  the  snow  cover.  Comparing  the  cases  with  and  without  snow 
cover  shows  that  the  low-loss  dry-snow  layer  can  enhsmce  the  backscattering  due  to  the 
introduction  of  more  scatterers  (ice  grains)  sind  a  better  impedance  match  between  the  air 
and  the  sea  ice  layer. 

The  above  results  have  shown  that  the  scattering  from  the  inhomogeneities  is  dom¬ 
inant  at  large  incident  angles  while  the  contribution  from  the  rough  surface  is  important 
at  small  incident  angles  for  co-polaxized  returns.  The  cross-polsirized  return  due  to  the 
inclusions  comes  from  the  first-order  term  under  the  distorted  Born  approximation  in  this 
model.  The  spheroidal  model  in  [2]  predicts  a  zero  Vcdue  of  cross-poleurized  backscattering 
coefficient  under  the  first-order  approximation  because  the  scatterers  are  vertically 
aligned.  In  the  next  subsections,  effects  of  correlation  lengths  and  polarization  signature 
of  sea  ice  are  considered. 


b.  Study  on  Correlation  Lengths 


SHAPE 

wm 

tz'  / iy' 

tz'ly'tz' 

Reference 

0.700 

0.250 

1.200 

2.800 

1.714 

4.800 

0.210 

Expanded 

1.050 

0.375 

1.800 

2.800 

1.714 

4.800 

0.709 

Contracted 

0.467 

0.167 

0.800 

2.800 

1.714 

4.800 

0.062 

Thickened 

0.467 

0.375 

1.200 

1.245 

2.570 

3.200 

0.210 

Thinned 

1.050 

0.167 

1.200 

6.287 

1.143 

7.186 

0.210 

Elongated 

0.572 

0.204 

1.800 

2.800 

3.147 

8.824 

0.210 

Shortened 

0.857 

0.306 

0.800 

2.800 

0.933 

2.614 

0.210 

Table  1.  Correlation  lengths  for  various  scatterer  shapes 

The  geometry  of  the  brine  inclusions  in  sea  ice  is  not  available  from  the  ice  charac¬ 
terization  data.  To  study  the  effect  of  the  scatterer  size  and  shape,  scattering  coefficients 
are  calculated  with  variable  correlation  lengths.  To  isolate  the  effect  of  correlation  lengths, 
bare  sea  ice  with  smooth  surface  configuration  in  Figure  4  is  considered  to  avoid  the  contri- 


15 


bution  from  the  snow  layer  and  rough  interfaces.  The  correlation  lengths  used  in  the  data 
matching  serve  as  a  reference.  The  scatterer  is  expanded,  contracted,  thickened,  thinned, 
elongated,  and  shortened  by  varying  the  correlation  lengths  as  shown  in  Table  1  while  the 
fractional  volume  is  kept  constant  in  all  cases. 

When  all  the  correlation  lengths  are  varied  by  the  same  factor,  the  scatterer  shape 
remains  unchanged  since  the  correlation  length  ratios  stay  the  same.  For  the  szune  frac¬ 
tional  volume,  the  effective  number  of  scatterers  decreases  as  the  scatterer  is  expanded  and 
increases  when  the  scatterer  is  contracted.  The  conventional  backscattering  coefficients  as 
well  as  the  correlation  coefficient  p  —  <rh.hvv/y/<^khO‘vv  [2]  between  the  horizontal  and  verti¬ 
cal  returns  are  calculated  for  the  volume  expansion  and  contraction  with  the  same  factor  of 
1.5^.  The  results  are  plotted  in  Figure  10.  For  <Thh  aiid  <7„,,  the  scattering  coefficients  are 
reduced  much  more  when  the  correlation  volume  of  the  scatterer  is  contracted  as  compzired 
to  the  increases  in  the  case  of  volume  expansion.  This  show  the  nonlineair  relation  between 
the  volume  scattering  and  the  scatterer  size.  In  this  case,  an  underestimate  in  the  scatterer 
volume,  compared  to  the  overestimate  with  the  szime  factor,  will  lead  to  more  variation  in 
the  co-polarized  backscattering  coefficients.  An  interesting  observation  is  the  insensitivity 
of  the  cross-polarized  ratio  e  =  ahvjo'hh  as  shown  in  Figure  10.  This  is  because  the  depo¬ 
larization  is  caused  by  the  non-spherical  shape  of  the  scatterer,  which  is  kept  unchamged. 
For  the  correlation  coefBcient  p,  the  results  show  that  the  horizontal  and  vertical  waves 
become  more  correlated  as  the  correlation  volume  increases.  For  the  larger  volume,  the 
sea  ice  medium  is  more  lossy  and  wave  path  is  effectively  reduced.  Consequently,  the 
anisotropic  effect  causing  the  decorrelation  between  the  horizontal  and  vertical  waves  is 
weakened. 

When  the  correlation  length  ratios  are  varied  while  the  scatterer  volume  and  the 
fractional  volume  are  kept,  the  scatterer  shaped  is  deformed  and  the  number  of  scatterers 
remains  unchanged.  For  the  same  vertical  correlation  length  £^t ,  the  horizontal  correlation 
lengths  are  now  modified  by  a  factor  of  1.5  to  thicken  and  thin  the  scatterer  as  indicated  in 
Table  1.  In  Figure  11,  the  co-polarized  returns  are  rather  insensitive  since  the  horizontal 
cross  section  of  the  scatterer  is  not  varied  =  0.175  mm^).  It  is  obvious  that  the 

cross-polarized  ratio  is  strongly  dependent  on  the  scatterer  shape.  The  depolarization 
effect  is  more  prominent  as  the  scatterer  is  further  deformed  from  the  spheroidal  shape. 
For  the  magnitude  of  p,  the  variation  is  larger  at  small  incident  angles  where  the  scatterer 
shapes  are  more  different  as  compared  to  the  cases  at  large  incident  angles.  The  phase  of  p 
started  from  0°  at  normal  due  to  the  azimuthal  symmetry  [20]  and  increases  with  incident 
angles.  Also,  the  absolute  value  of  the  phase  is  larger  for  the  thicker  scatterer  since  the 
anisotropy  is  stronger. 

While  the  correlation  length  ratio  lx>/ty>  and  the  fractional  volume  are  not  changed, 
the  scatterer  is  elongated  and  shortened  by  varying  the  vertical  correlation  length  by 
a  factor  of  1.5.  Figure  12  presents  the  results  for  these  cases  compared  with  the  reference. 
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There  are  large  changes  in  the  co-polarized  returns  due  to  the  corresponding  differences 
in  the  horizontal  cross  section  of  the  scatterer.  The  variation  in  the  cross-polarized  ratio 
is  small  since  the  horizontal  correlation  length  ratio  is  not  varied.  The  magnitude  of  p 
differs  more  at  large  incident  angle  where  the  shape  difference  is  more  significcint.  The 
phase  of  p  is  closer  to  zero  as  the  scatterer  becomes  shorter  and  the  sea  ice  medium 
approaches  the  isotropic  condition.  As  shown  in  this  subsection,  the  scatterer  volume  and 
the  horizontal  cross  section  affect  the  level  of  the  co-polarized  returns.  The  depolarization 
effect  characterized  by  the  cross-polarized  ratio  is  controlled  by  the  scatterer  shape.  The 
ellipsoidal  shape  of  the  scatterer  with  preferential  vertical  alignment  is  also  responsible  for 
the  anisotropy  of  the  medium  and  consequently  the  behavior  of  the  complex  correlation 
coefficient  p  between  the  horizontal  and  vertical  returns. 

c.  Pol2Lrimetric  Simulation 

To  investigate  how  polarization  signatures  relate  to  the  corresponding  covariance 
matrix  and  convey  information  regarding  different  sea  ice  configurations  emd  scattering 
mechanisms,  polaximetric  data  sure  simiilated  for  sea  ice  with  bare  surface  smd  with  snow 
cover  at  different  incident  angles.  The  configurations  under  consideration  are  in  Figure  4 
for  bare  sea  ice  with  smooth  surface,  in  Figure  6  for  bare  sea  ice  with  rough  surface,  and 
in  Figure  8  for  snow-covered  sea  ice.  Due  to  the  azimuthal  symmetry  of  the  scattering 
configuration,  the  covariance  matrix  has  the  form  of  [20] 

^  [•  1  0  py/j' 

C  =  <rhh  0  e  0  (32) 

.P*y/y  0  7  . 

where  the  co-polarized  ratio  is  7  =  <Tvvl<fhh'  The  elements  of  the  covariance  matrix  is 
used  to  calcidate  the  Mueller  matrix  to  obtain  the  signature  <r{ariPryOi,^i)  as  defined  in 
[2]  for  polarization  angles  and  a;  and  ellipticity  angles  fir  and  Pi  of  the  received  (r) 
and  the  incident  (t)  waves,  respectively.  For  a  co-polarized  signature,  =  oti  =  o:  and 
Pr  =  0i  =  P’ 

Normalized  co-polarized  signatures  denoted  by  <7n  for  bare  sea  ice  with  smooth  and 
rough  surfaces  at  incident  angles  of  20”  and  45”  are  presented  in  Figure  13.  At  ^ot  =  20”, 
the  signatures  for  smooth  sea  ice  due  to  volume  scattering  and  for  rough  sea  ice  due  to 
surface  scattering  are  similar  because  both  covariance  matrices  have  7  close  to  unity,  small 
e,  and  p  with  small  phase  sis  seen  in  Figure  14.  However,  the  magnitude  of  p  shown 
in  Figure  14  are  distinctively  smaller  for  the  volume  scattering  mechanism  with  smooth 
surface.  At  =  45”,  the  volume  scattering  is  dominant  in  both  smooth  and  rough  cases 
whose  correlation  coefficients  have  larger  phase  giving  rise  to  the  signatures  in  Figures  13c 
and  13d  with  more  distortion  compared  to  those  at  20”.  For  snow-covered  sea  ice.  Figure 
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15  reveal  that  |/9|  is  closer  to  unity  and  the  phase  is  smaller  at  incident  angles  of  20°  and 
45°.  Thus,  the  corresponding  signatures  for  the  snow-cover  sea  ice  in  Figure  15  have  much 
less  distortion  due  to  the  effects  of  the  rough  surfaces  and  the  isotropic  snow  layer  which 
mask  the  anisotropic  information  from  the  lower  sea  ice  layer  at  the  X-band  frequency. 
The  boundary  effects  due  to  snow  cover  can  also  be  observed  with  the  oscillation  in  the 
phase  of  p  in  Figure  14. 

5.  SUMMARY 

In  this  paper,  an  anisotropic  layer  model  with  ellipsoidal  scatterers  have  been  de¬ 
veloped.  The  top  layer  is  an  isotropic  medium  with  randomly  oriented  spheroids.  In  the 
anisotropic  layer,  the  ellipsoids  are  aligned  vertically  and  oriented  randomly  in  azimuthal 
directions.  Further  medium  complexity  such  as  multi-species,  size  and  shape  distributions 
have  not  been  incorporated  in  this  model  and  can  be  considered  for  future  model  develop¬ 
ment.  The  effective  peruxittivities  are  derived  with  the  extended  strong  fluctuation  theory 
and  the  polarimetric  sc  ttering  coefficients  are  calculated  under  the  distorted  Born  approx¬ 
imation.  In  this  irodel,  the  cross-polarized  return  is  obtained  from  the  first-order  distorted 
Born  approximation  due  to  non-spherical  shapes  of  the  scatterers;  however,  higher  order 
contributions  axe  ignored.  The  effects  of  rough  surfaces  are  estimated  with  the  incoher¬ 
ent  addition  method  and  higher-order  interactions  between  volume  and  surface  scattering 
mechanisms  are  not  included.  The  theoretical  results  compare  weU  with  the  measured  data 
for  thick  first-year  sea  ice  with  bare  surface  and  snow  cover.  The  effects  of  the  scatterer 
geometry  are  investigated  by  varying  correlation  lengths  which  change  the  shape  and  the 
volume  of  the  scatterer.  Polzirization  signatures  are  then  simulated  for  the  thick  first-yeeir 
with  and  without  snow  cover.  The  volume  and  surface  scattering  mechemisms  are  discussed 
in  terms  of  conventional  backscattering  coefficients  as  well  as  polarization  signatures  and 
correlation  coefficient  p  between  the  vertical  and  the  horizontal  returns. 
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Figure  Captions 


Figure  1.  Scattering  configuration  of  layer  media. 

Figure  2.  Illustration  of  a  sea-ice  horizontal  thin  section;  the  arrows  represent  random 
horizontal  orientations  of  crystallographic  c-axes  and  ellipsoidcd  brine  inclusions. 

Figure  3.  Geometry  of  an  ellipsoidal  scatterer. 

Figure  4.  Scattering  configuration  of  bare  thick  first-year  sea  ice  with  smooth  surface. 

Figure  5.  Backscattering  coefficients  from  bare  thick  first-year  sea  ice  with  smooth 
surface. 

Figure  6.  Scattering  configuration  of  bare  thick  first-year  sea  ice  with  rough  surface. 
Figure  7.  Backscattering  coefficients  from  bare  thick  first-year  sea  ice  with  rough  surface. 
Figure  8.  Scattering  configuration  of  thick  first-year  sea  ice  with  snow  cover. 

Figure  9.  Backscattering  coefficients  from  thick  first-year  sea  ice  with  snow  cover. 
Figure  10.  Effects  of  expanded  and  contracted  scatterer  volumes. 

Figure  11.  Effects  of  thickened  and  thinned  scatterer  shapes. 

Figure  12.  Effects  of  elongated  and  shortened  scatterer  shapes. 

Figure  13.  Polarization  signatures  of  bare  thick  first-year  sea  ice  with  smooth  and  rough 
surfaces  at  incident  angles  of  20°  and  45°. 

Figure  14.  Correlation  coefficient  p  of  thick  first-year  sea  ice  with  bare  smooth  surface, 
bare  rough  surface,  and  snow  cover. 

Figure  15.  Polarization  signatures  of  thick  first-year  sea  ice  with  snow  cover  at  incident 
angles  of  20°  and  45°. 
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ABSTRACT 

A  model  to  calculate  the  effective  permittivity  of  saline 
ice  under  thermal  variation  is  presented  in  this  paper.  The  model 
includes  multi-phase  inhomogeneities  with  multiple  species  char¬ 
acterized  by  orientation,  size,  and  shape  distributions.  The 
model  is  then  used  to  derive  the  effective  permittivity  as  a  func¬ 
tion  of  temperature  under  the  strong  fluctuation  theory  which 
is  e.xtended  to  account  for  the  complexity.  The  results  calcu¬ 
lated  from  the  model  are  compared  with  experimental  data  at 
4.8  GHz  for  saline  ice  grown  at  the  US  .Army  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory  (CRREL).  The  comparison 
between  measured  and  calculated  complex  permittivities  is  good 
for  the  imaginary  part  and  the  difference  is  within  10%  for  the 
real  part. 

INTRODUCTION 

Sea  ice  consists  of  crystalline  ice,  solid  salt,  liquid  brine 
inclusions,  and  air  bubbles.  Temperature  strongly  affects  sea 
ice  structure,  constituent  characteristics,  and  thus  its  electro¬ 
magnetic  properties.  The  constituent  phases  in  sea  ice  are  re¬ 
lated  thermodynamically.  Fractional  volumes  of  brine  inclusions 
and  air  bubbles  change  accordingly  as  the  temperature  of  sea 
ice  varies.  For  brine  inclusions,  the  size  distribution  has  been 
reported  to  follow  the  power  law  jlj  and  the  shape  has  been  ob¬ 
served  as  substantially  ellipsoidal  j2|  for  large  and  medium  sizes 
and  more  rounded  for  small  sizes. 

When  the  temperature  increases  during  a  warming  cy¬ 
cle,  the  shape  of  the  inclusions  becomes  less  ellipsoidal.  For  air 
bubbles,  the  shape  has  been  considered  as  rather  rounded  in  the 


form  of  spheroids  or  spheres,  which  does  not  change  as  much 
as  the  shape  of  brine  inclusions.  Non-spherical  inclusions  in  sea 
ice  can  have  preferential  alignment  in  the  vertical  direction  such 
as  in  columnar  ice  and  random  orientation  in  azimuthal  direc¬ 
tions.  In  this  model,  the  orientation  distribution  of  scatterers  b 
depicted  with  Eulerian  angles.  The  size  variation  of  scatterers 
in  a  species  is  described  in  terms  of  the  number  density  or  the 
fractional  volume  as  a  function  of  normalized  volumetric  sizes. 
For  various  shapes  of  a  scattering  constituent,  inclusions  with 
similar  shape  are  treated  collectively  as  a  subspecies.  Strong 
permittivity  fluctuation  theory  |3l  is  then  extended  to  derive  the 
anisotropic  permittivity  tensor  of  sea  ice  with  the  above  charac¬ 
teristics.  Theoretical  results  are  then  compared  with  measured 
data  for  saline  ice  at  different  temperatures. 

EFFECTIVE  PERMITTIVITY 

Consider  an  inhomogeneous  medium  consisting  of  vari¬ 
ous  scattering  species  embedded  in  a  background  medium.  First, 
a  species  is  defined  as  a  set  of  all  scatterers  with  the  same  permit¬ 
tivity.  Within  a  species,  scatterers  of  similar  shape  are  classified 
into  a  subspecies.  For  simplicity,  a  species  is  now  redefined  as  a 
set  of  all  scatterers  with  the  same  permittivity  and  shape.  Let 
ts  be  the  permittivity  of  the  background.  The  inhomogeneous 
medium  has  N  scatterer  species  of  ellipsoidal  shapes  (which  can 
be  reduced  to  spheroidal  or  spherical  shapes)  with  vertical  align¬ 
ment  and  random  azimuthal  orientation.  Species  i  has  permit' 
tivity  t,-,  and  a  fractional  volume  ft,  in  the  mixture.  An'isotropic 
effective  permittivity  tensor  of  the  inhomogeneous  medium  b  • 
sum  of  a  quasi-static  part  and  a  scattering-effect  part  in  the 
following  expression 
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where  auxiliary  permillivity  <j  and  dyadic  coefficient  S  are  de¬ 
rived  from  the  condition  of  secular-term  elimination,  and  the 
effective  dyadic  scatterer  is  approximated  as 


(W/]>m  =  H  /  dop.(a,/?,-y) 

£  {*0  [G,(r)i,,4.,«(r)  + 
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where  i  stands  for  species  i(s).  Quantity  u  is  the  normalized  vol¬ 
umetric  size  defined  as  the  ratio  of  the  volume  of  a  scatterer  over 
that  of  the  smallest  scatterer  of  the  same  permittivity  regardUss 
cj  shape.  The  integration  over  dv  accounts  for  the  size  distri¬ 
bution  in  the  size  range  i>,  of  the  scaitercrs  in  species  i.  The 
probability  density  function  of  orientation  is  p,(a,3,‘t]  where 
a,  fi,  and  •)  are  Eulerian  angles.  Gg  is  the  anisotropic  Green's 
function.  is  the  Fourier  transform  of  normalized  exponential 
correlation  function  with  correlation  lengths  f,i<,  and  <,i‘. 

is  local  variance  of  species  i.  The  condition  of  secular 
elimination  (j(r))  =  0  is  imposed  to  find 
and  then  the  integration  in  (2)  is  carried  out  to  derive  the  effec¬ 
tive  permittivity  tensor. 


PHYSICAL  DESCRIPTIONS 

At  microwave  frequencies,  saline  water  in  brine  inclu¬ 
sions  has  high  permittivity  which  strongly  affects  electromag¬ 
netic  properties  of  sea  ice.  Size  distribution  of  the  inclusions  has 
been  reported  to  follow  the  power  law  jlj.  The  power-law  distri¬ 
bution  in  terms  of  normalized  volumetric  size  is  described  by  a 
power-law  index  p  and  the  number  density  riQ  of  brine  inclusions 
of  smallest  size.  The  normalized  volumetric  size  is  defined  as  the 
ratio  of  volume  of  an  inclusion  over  the  volume  of  the  smallest 
one.  When  the  total  fractional  volume,  the  minimum,  average, 
and  maximum  sizes  of  the  brine  inclusions  are  measured,  p  and 
no  can  be  found  and  the  size  distribution  is  thus  determined.  .\t 
different  temperature,  volumes  of  individual  brine  pockets  vary. 

%  assuming  that  volumes  of  all  brine  inclusions  change  at  the 
same  rate  and  the  density  number  of  the  inclusions  is  conserved, 
the  power  index  can  be  shown  to  remain  unchanged  as  the  tem¬ 
perature  varies.  If  characterization  data  are  not  available  at  all 
'etnperatures  under  consideration,  the  index  obtained  at  a  tem¬ 
perature  may  be  used  at  a  different  temperature  when  there  is 
Severe  brine  loss  nor  deformation.  Otherwise,  size  measure¬ 
ments  are  necessary  to  determine  the  size  distribution.  For  air 
^  es,  the  size  distribution  is  derived  in  the  same  manner  or 
tained  from  section  images. 

.  l>e«n  observed  that  the  shape  of  brine  inclusions  is 

^ntially  ellipsoid  J2j.  From  ice'  thin  sections  i2',  only  inclu- 
small  sizes  are  seen  to  have  a  more  rounded  spheroidal 
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form.  Based  on  this  observation,  a  slowly  varying  logarithmic 
function  is  used  to  describe  the  shape  distribution  in  terms  of 
normalized  volumetric  size  and  axial  ratios  of  the  ellipsoids. 
When  the  temperature  is  increased  during  a  warming  cycle,  the 
shape  of  the  inclusions  becomes  more  rounded.  This  structural 
variation  in  sea  ice  is  characterized  with  reshaping  factors  op¬ 
erating  on  the  maximum  axial  ratios.  This  thermal  variation 
depicts  the  smoothing  process  which  makes  the  brine  inclusions 
approach  the  more  rounded  spheroidal  form  with  increased  tem¬ 
perature.  For  air  bubbles,  the  shapes  have  been  considered  as 
rather  rounded  in  the  form  of  spheroids  or  spheres.  Therefore, 
the  shapes  of  air  bubbles  do  not  change  as  much  as  those  of  brine 
inclusions  and  are  assumed  to  be  unaffected  by  temperature. 

Local  correlation  lengths  of  a  scatterer  can  also  be  es¬ 
timated  by  the  equivalence  condition  between  the  correlation 
volume  and  the  volume  of  the  scatterer  with  the  same  axial  ra¬ 
tios  l^j.  When  the  size  distribution  'is  specified,  local  correlation 
lengths  of  a  brine  inclusion  can  therefore  be  calculated.  Further¬ 
more,  correlation  lengths  also  vary  as  a  function  of  temperature 
since  the  size  of  an  individual  inclusion  changes  with  tempera¬ 
ture.  If  section  images  of  sea  ice  are  available  at  the  temperatures 
under  consideration,  sizes  of  scatterers  are  measured  to  obtain 
the  correlation  lengths.  If  the  size  measurements  are  only  made 
at  a  fix  temperature  T,-,,  correlation  lengths  at  a  different  tem¬ 
perature  T  can  be  estimated  for  undeformed  ice  with  low  brine 
loss.  Correlation  lengths  for  air  bubbles  can  also  be  obtained  in 
the  same  manner. 

The  constituents  in  sea  ice  are  ice,  solid  salt,  brine  inclu¬ 
sions,  and  air  bubbles.  For  the  background  ice,  the  real  part  of 
permittivity  is  not  sensitive  to  temperature;  however,  the  imagi¬ 
nary  part  is  dispersive  and  varies  with  temperature.  An  empiri¬ 
cal  formula  to  calculate  the  imaginary  part  of  ice  permittivity  as 
a  function  of  temperature  at  microwave  frequencies  is  given  in 
|5j.  For  brine  in  sea  ice,  empirical  formulas  to  compute  complex 
dielectric  constants  of  brine  in  terms  of  temperatures  and  mi¬ 
crowave  frequencies  are  reported  in  (6|.  To  determine  fractional 
volumes  of  brine  inclusions  and  air  bubbles  in  sea  ice,  Cox  and 
Weeks  i7|  have  provided  equations,  based  on  phase  equilibrium, 
in  terms  of  bulk  ice  density  and  salinity  for  the  temperature 
range  of  2  C  lo  -30  C.  The  fractional  volume  of  solid  salt  is 
low  and  will  be  neglected  in  the  calculation.  From  characteriza¬ 
tion  measurements,  input  parameters  to  the  sea  ice  model  can 
be  obtained  as  discussed  in  this  section. 

RESULTS  AND  DISCUSSIONS 

Relative  permittivities  of  ice  were  measured  ;2i  at  CR- 
REL  by  using  the  transmission  method  with  wave  normally  inci¬ 
dent  on  ice  slabs;  thus,  only  t,//,  corresponding  to  the  ordinary 
wave  were  obtained.  Ice  slab  85-3  considered  here  was  taken 
from  a  saline  ice  sheet  to  a  laboratory  environment.  Measure¬ 
ments  of  permittivity  and  physical  parameters  were  made  during 
progressive  warming  of  the  ice  slab  from  -32'’C  to  -2°C  at  4.8 
GHz.  During  the  process,  there  w'as  some  brine  loss  which  can 


also  be  incorporated  in  the  calculations.  Results  are  shown  in 
the  enclosed  figure  for  the  imaginary  part  of  the  relative  permit¬ 
tivity  of  the  saline  ice.  When  the  temperature  is  increased,  the 
smoothing  effect  of  rounded  inclusions  is  in  competition  with  the 
increasing  effect  of  higher  fractional  volumes  to  render  the  slowly 
increasing  trend  in  the  imaginary  part.  The  drop  at  tempera¬ 
tures  higher  than  -8'C  can  be  explained  by  the  brine  loss  and 
the  kink  observed  at  about  -22'C  is  due  to  the  phase  change  as¬ 
sociated  with  the  precipitation  of  sodium  chloride  dihydride.  As 
seen  from  the  figure,  the  theoretical  curve  explains  all  the  trends 
observed  in  the  experimental  data  for  the  imaginary  part  of  the 
relative  effective  permittivity.  The  real  parts  of  the  permittivity 
is  rather  monoton  ically  increased  from  low  to  high  temperatures. 
The  calculated  values  are  lower  but  are  within  1055  of  the  mea¬ 
sured  data. 


Temperolure  (®C) 


SUMMARY 

A  model  which  accounts  for  the  complexity  of  the  multi¬ 
phase  inhomogeneities  with  multiple  species  characterized  by  ori¬ 
entation,  size,  and  shapes  distributions  has  been  presented  in 
this  paper.  Anisotropic  effective  permittivity  of  the  inhomoge¬ 
neous  medium  is  derived  with  the  extended  strong  fluctuation 
theory.  The  thermodynamics  of  the  constituents  in  sea  ice  is 
also  included  in  the  model  to  incorporate  the  thermal  variations 
in  the  ice  structure  and  properties.  A  comparison  of  theoretical 
results  with  measured  data  for  effective  permittivities  of  saline 
ice  is  made  and  the  trends  observed  from  the  experimental  mea¬ 
surements  are  explained. 
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Abstract 

Earth  tenains  are  modeled  by  a  two-layer  configuration  with  a  random  permittivity  described 
by  a  three-dimensional  correlation  function  with  horizontal  and  vertical  correlation  lengths  and 
variances.  Using  the  wave  theory  with  Bom  approximations  carried  to  the  second  order,  this  model  is 
applied  to  derive  the  polarimetric  backscattering  coefficients  of  the  Mueller  and  covariance  matrices. 
Fl'om  a  physical  point  of  view,  the  Bora  first-  and  second-order  approximations  account  for  the  single 
and  double  scattering  processes,  respectively. 

For  the  isotropic  random  permittivity  configuration,  five  out  of  the  nine  elements  of  the 
covariance  matrix  are  zero  under  the  Born  first-order  approximation.  For  the  uniaxial  tilted  random 
permittirity  case,  the  covariance  matrix  does  not  contain  any  zero  elements.  To  account  for  the 
azimuthal  randomness  in  the  growth  direction  of  leaves  in  tree  and  grass  fields,  an  averaging  scheme  is 
developed  in  which  the  backscattering  coefficients  are  averaged  over  the  azimuthal  direction.  In  this 
case,  the  covariance  matrix  is  described  by  four  zero  elements  though  the  tilt  angle  is  different  from 
zero.  Applying  the  Born  approximation  to  the  second-order,  the  covariance  matrix  is  computed  for 
the  isotropic  and  the  uniaxial  untiited  random  permittivity  configurations.  The  covariance  matrix 
has  four  zero  elements,  and  a  depolarization  factor  is  obtained,  even  for  the  isotropic  case. 

To  describe  the  effect  of  the  random  medium  on  electromagnetic  waves,  the  strong  permittiv¬ 
ity  fluctuation  theory,  which  accounts  for  the  losses  due  to  both  of  the  absorption  and  the  scattering, 
is  used  to  compute  the  effective  permittivity  of  the  medium.  For  a  mixture  of  two  components,  only 
the  frequency,  the  correlation  lengths,  the  fractional  volume,  and  the  permittivities  of  the  two  con¬ 
stituents  are  needed  to  obtain  the  polarimetric  backscattering  coefficients.  Theoretical  predictions 
are  illustrated  by  comparing  the  results  with  experimental  data  for  vegetation  fields  and  sea  ke. 
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Introduction 

For  many  types  of  earth  terrain,  the  scattering  effects  due  to  medium  inhomo- 
geneities  play  an  important  role  in  the  determination  of  radar  backscattering  coeffi¬ 
cients.  The  volume  scattering  properties  of  a  medium  have  been  successfully  modeled 
by  the  random  permittitity  fluctuations  [1].  The  random  medium  is  characterized  by 
a  background  permittivity  and  a  random  permittivity  fluctuation  modeled  by  a  three- 
dimensional  correlation  function  with  horizontal  and  vertical  correlation  lengths  and 
variances. 

The  objective  of  this  paper  is  to  develop  a  mathematically  rigorous,  fuUy  po- 
larimetric  model  to  compute  the  covariance  matrices  in  the  backscattering  direction 
for  different  kinds  of  earth  terrains  with  a  two-layer  configuration  medium  with  ei¬ 
ther  an  isotropic  or  an  anisotropic  random  permittivity.  Of  all  the  analytical  wave 
approaches  to  the  study  of  scattering  by  random  medium,  the  Born  approximation 
is  the  simplest.  An  integral  equation  is  first  formed  for  the  electric  field  using  the 
unperturbed  Green's  function  in  the  absence  of  permittivity  fluctuations.  The  inte¬ 
gral  equation  is  then  solved  by  iteration  with  the  iteration  series  known  as  the  Born 
series.  With  both  the  first  and  second  terms  of  this  series  which  respectively  describes 
a  single  and  a  double  scattering  processes,  the  polarimetric  backscattering  coefficients 
are  computed. 


Scattering,  Mueller,  and  Covariance  Matrices 

Consider  an  electromagnetic  plane  wave  propagating  in  the  direction  k. 

£  =  (£»h -b  £.i;)  e'*  '  (1) 

where  is  the  horizontal  component,  £,  is  the  vertical  component,  h  and  v  are  two 
orthogonal  unit  vectors  satisfying  the  relation  v  x  h  =  k.  Consider  an  incident  plane 
wave impinging  on  a  particle.  This  particle  will  radiate  and  the  scattered  electric 
field  E,  can  be  related  to  the  incident  electric  field  £,•  by  the  scattering  matrix  (2): 

(!::)=(?? 

The  elements  of  the  scattering  matrix  are  complex  and  the  relation  V'/f  =  H\'  holds 
[1]  in  the  backscattering  direction  and  for  reciprocal  media. 

Instead  of  relating  the  scattered  electric  field  to  the  incident  electric  field,  the 
modified  Stokes  vector,  defined  by  (l),  cau  also  be  used  to  relate  the  scattered  Stokes 
vector  to  the  incident  Stokes  vector: 

=  /v  =  ^  b’ =-£«(£»£;)  V  =  -/m(£*£;)  (3) 

V  ^  V  V 


This  definition  is  chosen  since  the  first  two  modified  Stokes  parameters  represent  the 
intensity  of  the  wave  In  the  horizontal  and  vertical  polarizations: 


(4) 


In  terms  of  the  elements  of  the  scattering  matrix  (2),  the  Mueller  matrix  M  relates 
the  scattered  to  the  incident  Stokes  vectors  according  to 

7.=^-  7.  (5) 


is  found  to  be 

-tmiHH-HV)  \ 

\VHf  \VV\^  Re{VHV\'*)  -Im(VH-VV)  \ 

iRt{HHVH')  2Rt{BVVV)  Re(Utt  VV*  +  HV -VH')  VV*  -  BV  Vff*)  j 

7tm{BB-VH*)  2lm(BVVV)  Im{BB  ■  VV  +  BV  VB*)  Re{BB  VV’  -  BV  VB*)  J 

(6) 

where  all  the  Mueller  matrix  elements  are  real  numbers.  From  the  four  ele¬ 
ments  of  the  scattering  matrix  (2),  the  covariance  matrix  may  also  be  defined.  In 
general  the  co\'ariance  is  a  complex  matrix  of  sixteen  elements.  In  the  special  case 
of  the  backscattering  direction  and  for  a  reciprocal  medium,  the  covariance  matrix  is 
transformed  to  a  matrix  of  nine  elements  as  follows 

^  (HH\  {  HHHV  HHVV\ 

C=^\hv\  {HH"  HV  VV  )  =  I  ffV  .  HH*  \HVi^  HV-VV 
\VV  I  \VV.HH*  VVHV  \VV\^  ) 

(7) 

By  comparing  the  scattering  (2),  the  Mueller  (6),  and  the  covariance  matrices  (7), 
it  is  easily  seen  that  all  the  polarimetric  information  may  be  expressed  in  any  of 
these  matrices.  The  covariance  and  Mueller  matrices  can  be  transformed  to  each 
other  since  they  contain  the  same  information.  However,  from  these  tn-o  matrices, 
the  scattering  matrix  cannot  be  fully  retrieved  because  the  absolute  phase  factor  is 
not  taken  into  account  in  either  the  covariance  matrix  or  in  the  Mueller  matrix. 


In  our  study  of  microwave  remote  sensing,  the  random  medium  model  is  used 
to  simulate  the  permittmty  of  the  earth  terrain.  This  implies  that  the  simple  case 
of  a  single  panicle  can  no  longer  be  considered  but  the  ensemble  average  of  the 
electric  field  needs  to  be  calculated.  The  polarimetric  bistatic  scattering  coefficient 
^*j)  is  defined  by  |3j 


^*j)  —  lim 


4xt* 

Acos9oi 


<  E^E^,  > 

EriEl 


lim 


47rr» 

A  COS  6tn 


<  flT  •  VK*  > 


(8) 
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where  A  is  the  area  illuminated  by  the  radar,  r  is  the  distance  between  the  radar 
and  the  observed  terrain,  and  $oi  is  the  polar  incident  angle.  Furthermore,  it,-  and 
kt  are  the  incident  and  scattered  wave  vectors,  respectively.  The  subscripts  s  and 
i  stand  for  scattered  and  incident  waves,  respectively.  Finally,  fir  and  i//c  are  the 
four  components  of  the  scattering  matrix  where  ft,  r,  v,  and  k  represent  either  the 
horizontal  h  or  vertical  v  polarizations. 

The  polarimetric  Wkscattering  coefficient  is  defined  j3]  from  the  bistatic 
coefficient  when  k,  =  -ki 


—  ki}  COsBoi  — 


<  Eu 


EhE. 


El,  >  4Tr* 

=  lim 


•*—  A 


<  flT  VK  > 


(9) 

The  polarimetric  backseat tering  coefficients  are  complex  except  the  three  “tra¬ 

ditional”  real-\’alued  backseat  tering  coefficients 


<r.,«  =»  <Txx  (10) 


Similarly,  the  covariance  matrix  (7)  constituted  by  the  different  polarimetric  backseat- 
tering  coefficients  can  be  derived 


C 


^  Ohh 

Ohkhx 

EV 

■  ( HH*  HV  VV  )  = 

<hhx 

OKx 

1 

[vvj 

\olkxx 

<^kxxx 

/ 

(11.) 


Therefore,  the  complete  polarimetric  information  can  be  obtained  by  computing  only 
six  backscatiering  coefficients  in  which  three  are  real  (tr**,  a*,,  and  <r„)  and  three 
are  complex  (ohhhx,  ffkhTx,  Okxxx)- 


Formulation 

Consider  an  electromagnetic  plane  wave  impinging  upon  a  two-layer  random 
medium  configuration  as  shown  in  Fig.  1.  The  top  laj-er  represents  free  space  with 
a  permittivity  cq  while  the  bottom  layer  has  a  deterministic  permittivity  t2-  In 
between,  a  layer  with  thickness  d  and  extending  infinitely  in  the  lateral  direction  is 
characterized  by  a  random  permittivity  <i(f) 


«i(r)  =  fi,„  +  f  j/(r)  (12) 

This  random  variable  has  a  mean  value  <i,„  and  a  small  fluctuating  term  fi/(r)  whose 
average  is  zero.  The  wave  equations  for  both  regions  may  be  written  as 


V  X  V  X  Eo(f)  -  koEo(r)  =  0 


(13) 
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V  X  V  X  £,(f)  -  =  Q(V^i(r)  (14) 

vkbtre  £o(f)  and  £i(f;  denote  the  electric  fields  fluctuation  in  region  0  and  region  1, 
respectively. 


z 


Layer  With  Random  Permittivity 


z  >*d 


Figure  1:  Isotropic  two-layer  configuration. 


Funhermore,  the  following  variables  have  been  defined 

A-J  =  =  w  Vim  <?(?)  =  wVi/(’‘)  (15) 

With  the  effective  source  term  Q(f)£i(r),  (13)  and  (14)  can  be  expressed  in  integral 
form  as  follows 

^o(?)  =  £^(5=)  +  / d’ri  Goi(?,?i)  Q(ri)  Ei(fi)  (16) 

£,(f)  =  £;®’(f)^  |rf’r,  S„(f.  ?,)(?(?,)  •£,(?,)  (17) 

where  £j***(r)  and  are  the  specular  or  coherent  electric_fields  in  the  absence  of 

random  permittivity  in  region  1.  Furthermore  Goi(r,fi)  and  Gii(?,fi)  are  the  dj'adic 
Green’s  functions  for  observers  in  regions  0  and  1,  respectively,  with  the  source  in 
region  1.  Solutions  to  these  two  coupled  equations  can  be  obtained  by  iteration. 
Substituting  (17)  into  (16),  the  total  electric  field  in  region  0  is  found  to  be  in  form 
of  the  Neumann  series 

■>I 


(18) 
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where  the  n‘*  order  field  is  given  by 

=  Jd^fi'--d^f,Go\{r,fi)Q(fi]-Gi\(fi,fi]Q{fj)“-G\t{fn-i,fn)Q(f»)‘£\  (^«)  (19) 

The  so-called  “Born  first-order  approximation”  is  obtained  by  letting  n  =  1  in 
the  previous  equation.  Hence 

(20) 

and  the  “Born  second-order  approximation”  is  derived  by  letting  n  =  2 

d%  Soi(r,J=i)  (?(f,)  Sn(ri.f2)  (21) 


Bora  First-Order  Approximation 
a)  Isotropic  random  medium  model 


Consider  an  electromagnetic  plane  ^•ave  impinging  upon  a  two-layer  random 
medium  configuration  as  shown  in  Fig.  1.  For  an  incident  electric  field  in  region  0  of 
the  following  form 

So.(r)  =  +  £‘or.r(-ito=.)c-*-*]  (22) 

where  h{-ko,i)  and  v{-kQ,i)  are  respectively  the  horizontal  and  vertical  unit  vectors 
in  region  0  associated  with  down-going  waves,  the  coherent  electric  field  in  the  random 
medium  is  obtained  (4] 


£1“^?)  =  {£b*.  -h 

-I-  £o..-  [A..ti(l:u.)e‘*>“‘  + 
in  which  the  following  coefficients  are  defined  [3] 

ko  loii 


>-e 


Aa.  =  ^  £i2. 

„  _x„, 

D,i  =  1  +  RouRiTi  e*’**"" 
p  — 

Xabi  =  1  +  RcU 


A„  =  f 

M  CJi 
kt  You 
“  ki  Fx 

^2.  =  l  +  So,.Sij.e*“-‘' 

SCfrlkaxt  “  ^ak^xi 
«*•  ~  — ; — ; — 

Yabi  =  1  +  Sabi 


(23) 


(24) 


As  can  be  seen  in  (23),  the  A  coefficients  are  associated  with  the  up-going  waves  while 
the  B  coefficients  describe  the  down-going  wa\es.  Under  the  far-field  approximation, 
the  Green's  function  in  (20)  takes  the  following  form  in  the  backscattering  direction 

=  . 

Goi(r,f,)  =  A(-l:o„)Soih(f,ri)  -r  f(-Ioxi)G|)ir(r,f,)  (25) 


4;rr  t  J 

Gow(F.?j)  =  ^  [A..i'(fcw. )€•*>'-••  -}■  B,iii-kxxi)e-'^'"‘'] 


(26) 

(27) 


For  each  polarization,  both  Goi(r,f|)  and  £,  (f)  are  constituted  by  a  down-  and  an 
uphgoing  waves.  Therefore,  £^**(r)  is  formed  by  the  sum  of  four  different  terms  as 
shown  in  Fig.  2.  All  the  multiple  rcBections  occurring  at  the  boundaries  are  also 
incorporated  in  the  model  and  the  backscattered  wave  is  due  to  the  single  scattering 
process  under  the  Bom  first-order  approximation. 


r 


Figure  2:  Physical  interpretation 


bVfk  wviu 


•  •• 


(’)■  componem 

inriHpnt  \  •  ^  Deeded.  Furthermore,  the  polarization  of  the 

the  FounTtrfr^fl"  Vf  H  ^onzontaJ  or  vertical.  With 

statistical  homogeneTus  raldU“medTum*‘‘°“  P^mittmty  fluctuation  for  a 


<  Q(^i)Q(r2r  >= 


(28) 
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the  ensemble  average  for  the  first-order  backscattered  intensity  takes  the  following 
form 

<  4’J(r)Eil>(fr  >  =  / d^p  «(^)e-^ 

[Go,^(r.ri) .  [5o:.(f,?2)  •  (29) 

where  t  and  k  represent  either  the  horizontal  or  vertical  component  of  the  transmitted 
w’ave  and  fi  and  v  represent  either  the  horizontal  or  vertical  component  of  the  r.xeived 
wave.  FVom  (23),  the  following  term  has  been  defined 

^il’(?)  =  E  Q(pA:i„)e'>*-'e'V-?  (30) 

p=-i 


where  a  stands  either  for  k  or  v,  and  p  describes  either  up  or  down-going  waves.  In 
a  similar  manner, 

Soi„(f,r,)  =  ^  A«p.  a(p^-l,i)e‘'*^'•^‘c*‘''■^‘  (31) 

P=-i 

where  the  following  definitions  have  been  used 


Xku  —  Aki  Xk-u  =  Bki  X,u  =  A„i  Xv-u  =  (32) 

Substituting  (30)  and  (31)  into  (29),  a  sum  of  16  terms  is  obtained 

<  £op  (r)£jL‘(f)*  E  A(pfci«)  •  HM] 

(X.wA'„.  Hrki,i)-k{skui)Y  ]  d'f.d'fj  j  i{P) 

As  can  be  seen,  all  the  polarization  information  is  contained  outside  the  integral. 
Also,  the  intensity  is  expressed  as  a  sum  of  sixteen  terms  since  the  electric  field  is 
composed  of  four  elements.  With  the  two  foLowing  scalar  quantities 

=  •Ypp.'XVYVX;.,  (34) 

■  r(9i-u-.)j  HrA-,..,)  ■  k\skui)T  (3ol 

(33)  can  be  written  as 


<  4i’(f)4:’(r)-  >  = 


16jr*r* 


^2  4’;r,nr~. 


J d^p  *(5)  [e’^* 


.+«il4.)i|g2ii^P,j 


(36) 
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The  integrals  in  (36)  axe  solved  by  decomposing  the  variables  of  integration  into  their 
lateral  and  vertical  components.  By  using  the  property  of  the  Delta  function  and 
introducing  the  illuminated  area  A,  the  integration  over  the  lateral  variables  p,,  pj, 
and  yields 


<  >  = 


4r2  ' 


In  order  to  solve  the  last  set  of  equations,  let  us  define  first  the  two  complex  quantities 


a  =  pki,i  +  qki,i  b  =  (rfci„-  +  sA'i.ij*  (38) 

Hence,  the  integral  part  of  (37)  may  be  written 

J  dz,dz2  I  dp.  $(^^  =  2it^-,/3,)  (39) 

Due  to  Cauchy’s  theorem,  the  integral  over  P.  yields  the  foiloTang  expression  by 
distinguishing  two  cases  depending  upon  the  sign  of  zi  -  Z2 


J dP.  i{-2k..i,P.) 


2jr:  Res^2kpi,0*)e~^^*'~’^^ 
-2Tti  Res^(2kpi,P~)e~'^‘^’'~‘^^ 


if  ii  <  Z2 
if  2i  >  Z2 


(40) 


where  it  b  assumed  that  the  spectrum  density  ^(p^,p.)  has  only  one  pair  of  complex 
conjugate  poles  at  P.  =  pf,  and  i(p.,.P,)  vanishes  everywhere  on  the  infinite  circle. 
For  the  random  medium  of  infinite  lateral  dimension  and  vertical  extension  from 
-  =  -  d  to  r  =  0,  (39)  becomes 


2x1  Res^{p^  =  2k^,  p* )  dzi  dr,  e“' ) 

-  2xt  Res^'P.,  =  2fp.. P;)  dzi  dzj  )  (41) 

For  Res{p.,  =  2kpi,P~)  =  -Res0.,  =  2kp{,p^),  these  integrals  may  be  easily  solved 
to  get 


2xt  Res^p^  =  ‘^k..i,p*) 


1  _  g-(«-^7)<t 

(a-Pt)(b-pt) 


I  _ 

{a-b)(b-P:) 


1  -  c’(*~^‘')‘' 

'^{a-b){b-P:)  ~  {a- P7){b-p7)  (a-p.-)(b-p-) 


(42) 
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To  simplify  the  expression,  let  us  define  the  following  term 
1  _  g— (“-A*  y 


_ 


1  _  1  _  g-«(a-6)<i 

+ 


(a  -  0t)  (*  -  0t\  {<^  -  i)  (*  -  Pt)  (<i-b){b-  0:) 


{a-0:){h-0:)  {a-0:)(b-0:) 


(43) 


FVom  (38),  a  is  different  from  6  for  most  cases  because  of  the  complex  conjugation. 
However,  for  the  special  case  when  a  =  b  —  0,  (43)  may  be  written  as 


1  -  id  id  1  - 

Pf  ^  Pt  ‘  P:  ^  P:^ 


(44) 


Substituting  (43)  and  (44)  into  (37)  and  using  (9),  the  polarimetric  backscattering 
coefficients  are  found  to  be 


o,,..  =  =  2k^,0:)  Y.  («) 


Only  three  polarimetric  backscattering  coefficients  are  different  from  zero.  This  comes 
from  the  factor  in  (35).  They  are  zero  due  to  the  definition  of  the  unit  vectors 
{3j.  Hence,  the  final  results  are 


OKK  =  iResH0^  =  2k^,  0^ )  £  ^Zhh 

<T„ = iResi(0,  =  2k^. $:)  Y 

- 1 

c,,„^27:Hu*fP\€,„\UResi(3^^2k^,0:)  Y  k::,  m 


^Ar  *  ^AtAt  *  ^AAAt  •”  ^hvvv  —  0  (4  0 

For  an  autocorrelation  function  for  the  random  permittivity  (12)  of  the  form 


i«i  -»a'  ln-«.  Ui-i-' 

<  <i/(ri)<i/(r2)'  >=  flfiml  e  c  ■'  e"  '* 


(48) 


where  f  is  the  variance.  Ip  and  f,  are  the  lateral  and  vertical  correlation  lengths, 
respectively,  the  spectrum  density  (28)  is  given  by 


^P.,Py,P.)  = 


V** 


(l  +  /?Ji2)  (1  +  ^/2) 


(49) 
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The  spectral  density  goes  to  zero  as  ^  tends  to  infinity  satis^-ing  the  requirement  lor 
the  complex  integration  (40).  The  poles  0,  of  the  spectrum  density  are  obtained 


.50) 


Furthermore,  the  associated  residues  are  found  to  be 


Resi(0„0y,0f) 


'2i:H  {1+0111)  (l  +  i^a) 


b)  Anisotropic  random  medium  model 

In  order  to  describe  earth  terrains  exhibiting  an  anisotropic  behavior,  the  ran* 
dom  permittivity  is  chosen  to  be  uniaxiai  in  a  tilted  coordinate  system  |6j.  Earth 
terrains  such  as  sea  ice  and  com  fields  display  such  a  property.  From  (12),  the  random 
permittivity  is  writteu  in  a  tensor  form 

?i(f‘)  =  ?ifB +  ?!/(»*)  5®2) 

where  ti/(f)  is  a  fimction  of  position  characterizing  the  randomly  fluctuating  com¬ 
ponent  of  the  permittivity.  For  the  random  medium,  it  is  assumed  that  both  fi„ 
and  fi/(f}  are  uniaxial  in  a  coordinate  system  xr/z'  tilted  respectively  at  an  angle  V. 
xl-’f  around  the  r  axis  as  represented  in  Fig.  3.  In  this  tilted  coordinate  st-stem.  the 
permittivity  tensor  «i„  takes  the  following  form 


In  the  coordinate  system  xyz,  the  permittivity  tensor  can  be  expressed  as 

/«u  0  0  \ 

tlm  ~  I  ®  fZZ  <23  1  (^) 

\  0  <32  <33/ 

where 


<11  =  <i 

<22  =  <1  cos*  Ip  +  sin*  V 
<23  =  <32  =  (<u  -  <i)sint(>cosv'’ 
<33  =  <1  sin*  ip  4-  cos*  V' 


(55) 


7 


Figure  3:  Geometrical  configuration  of  the  pennitthity  tensor  in  an  anisotropic 
random  medium. 


Similar  equations  nay  be  written  for  fi/(r)  with  the  tUt  angle  ti’/-  From  (20),  the 
“Bom  first-order  approximation”  yields  for  the  anisotropic  case 

=  /  d’fi  Soi(F,  r, . .  5(Fi)  •  (56) 

where  Q(f)  =  u;*/i?i/(r).  Using  Euler  notation  and  decomposing  the  electric  field  in 
horizontal  and  vertical  components,  the  ensemble  av-erage  for  the  first-order  backscat- 
tered  intensity  may  be  written  as 

<  4i'(r)£ov  (r)’ >  =  Y,  f  [Goi„(f,F,)E}“i(r,)] 

{Coiwfr.rilE',®;,!?,)]  <  Q;*(fj)Qi*„(r2)  >  (oTi 

where  the  average  of  the  components  of  the  electric  field  has  been  computed  by  defin¬ 
ing  ft  and  V  as  the  received  polarizations  and  r  and  k  as  the  transmitted  polarizations. 

As  sho^-n  in  Fig.  4,  a  horizontally  or  a  vertically  polarized  incident  electric  field 
u-ill  generate  both  an  ordinary  and  an  extraordinary  waves,  due  to  the  anisotropic 
permittivity  in  region  1.  defined  by  their  respective  waw  numbers  |6] 


(58) 


Figure  4:  Ordinary  and  extraordinary  electromagnetic  waves  excited  by  the  ani¬ 
sotropic  random  medium,  (a)  Horizontally  polarized  incident  electric 
field,  (b)  Vertically  polarized  electric  field. 
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The  components  of  the  unperturbed  electric  field  ^***(r)  and  the  Green’s  function 
Goi(r,ri)  may  be  written  as  (7) 

i:  (59) 


Goio(?.fi)  =  ^  E  (60) 

p^o.-<neu^ed 

where  a  stands  for  either  h  or  v,  and  p  describes  either  an  ordinary  up-going  wave  (o), 
an  ordinary  down-going  wave  (-o).  an  extraordinary  up-going  (eu),  or  an  extraordi¬ 
nary  down-going  wave  (ed).  Both  the  incident  electric  field  and  the  Green’s  function 
are  formed  each  by  four  terms  meaning  that  backscattered  electric  field  computed 
with  Born  first-order  approximation  will  be  constituted  by  sixteen  terms.  Further¬ 
more,  the  coefficients  X^pi,  the  unit  vectors  p{k^,i),  and  the  wrave  numbers  itf„-  take 
on  the  following  expressions  depending  on  the  value  of  p 


P  =  O  =»  *rz.  =  kui 

Xapi  —  Aooi 

p  =  ~o  ^  Irf,,-  =  - k^^- 

Pik^x^)  =  K-ktJ 

Xapi  =  Baai 

p  =  «u  =»  =  kl,i 

P(*?ri)  =  i(iS.) 

Xapi  =  Aa« 

p  ~  ed^  k^^  —  kl^ 

Xapi  — 

(61) 

Hence,  substituting  (39)  and  (60)  into  (57)  yields 


E 


<  4i’(?)4l'(r)’  >=  £o..£o‘«  E 

[Xpp^Xr^  p(kiJi .  [X^X^  HkU,  •  HkLUy  j 


(62) 


As  can  be  seen,  all  the  polarization  information  is  contained  outside  the  integral. 
The  parameters  p.  q.  r,  and  s  can  take  the  values  o,  -  o,  eu,  or  ed  as  described  in 
(61).  The  intensity  is  hence  given  b>'  a  sum  of  236  terms.  For  each  of  these  terms,  81 
products  (sum  over  jklm)  have  to  be  computed  due  to  the  anisotropic  beharior  of 
the  random  medium.  In  order  to  simplify  this  result,  let  us  define  the  two  following 
scalar  quantities 

KrS,  =  (63) 


=  ip(i-f„);  •  ^{kiM  if(i:i„)/  •  Hkuuv 


(64) 
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Hence,  (62)  may  be  written 


(65) 

By  following  a  similar  procedure  as  in  the  Isotropic  case  described  in  (37)  through 
(45),  the  parameters  a  and  6  are  first  defined  as 


(66) 


By  assuming  a  spectrum  density  with  n  pairs  of  complex  conjugate  poles 

at  0,  =  0f^  and  defining  the  following  factors 


/rr'  = 


1  _  ^-«(«-tf«)4  j  _  g-i(a-6)4 


(67) 


ic-0rJ^b.0tJ~  {a-.  b)(b-0tJ 

1  -  e"*'*“*^  «‘(*“^*'«)4 

(a -6)  (6-/?-)  ■  [a-0:„){b-0:J  (a  -  ,3- )(6  -  .i?- ) 
the  six  polarimetric  backscattering  coefficients  (11)  are  finally  found  to  be  (6] 


/nr‘  = 


=  27tVVi<ii’  x:  4-”;;,  S  njIJA 
iX  {fies$>«„(3,  =  2k^,0:j  -  Res^iUmid,  =  )  CT*} 


(69) 


Compared  to  the  results  of  the  isotropic  case  shown  in  (46)  and  (47),  there  is  no  zero 
elements.  This  is  due  to  the  anisotropic  behavior  and  the  excitation  of  both  ordinary 
and  extraordinary  waves. 


Born  Second-Order  Approximation 


For  this  case,  only  the  isotropic  random  pennitti\ity  will  be  considered.  FVom 
(21),  the  second-order  scattered  electric  field  is  given 

^0  *(^)  =  J  Goi{f,ri)  Q(fi)-Gii(fi,f2)  Q{f3)-E^i\f2)  (70) 


where  £,_(r)  and  Goj(r.fi)  are  given  by  (30)  and  (31),  respectively.  The  Green’s 
function  Gii(fi,f2)  relates  the  field  observed  in  region  1  at  a  position  fj  due  to  a 
source  in  region  1  located  at  fj.  It  takes  the  following  form  from  [9] 


where 


-2) 

ln(Sp.  21.^2) 

1  f  ■*/  (?!-?,)  f  -1  >  ^2 

J  1  2l  <  r2 

(71) 

Z  C^(^e)6(«ki.)a(-rAu) 

(72) 

Z  o{9fcir)a(-rfcu) 

(73) 

Here  q  denotes  the  polarization  direction  along  A  or  r,  and  the  coefficients  C  take 
the  following  values  [3] 


,(r^)  = 

C*^n(fcp)  =  Ctx-xiM,)  = 

Cr-n(fcp)  =  C<  .,(fcp)  =  SioSirf’^'-VF* 


Clx-xiK)  =  CUxx^,)  =  1/Z)2 
CU-xKK)  =  C^,.,(fcp)  = 

C,-,(X>)  =  C<.,.,(Iv)  =  S,o/f2 

(74) 


The  physical  meaning  of  the  Born  second-order  approximation  (70)  are  interpreted 
based  on  the  fact  that  Sod^-ri).  Gn(Fj.r2),  and  are  composed  of  a  down- 

and  an  up>-going  wave  for  each  polarization.  Therefore.  £”’(?)  is  formed  by  the  sum 
of  sixteen  different  terms  as  shown  in  Fig.  o.  All  the  multiple  reflections  occurring 
at  the  boundaries  are  accounted  for  and  the  backscattered  wave  is  due  to  the  double 
scattering  process  under  the  Bom  second-order  approximation. 

With  (30),  (31),  and  (70),  the  ensemble  average  for  the  second-order  intensity 
is  given 

<  4^’(7)  £5!’(?)*  >=  /d’F.d’rjd’rjd’f,  <  (?(?i)G{f2)Q(f'3)’G(r4)*  > 

[So,p(f,r,) .  ^„(f,,f2)  •  rl®;(f2)]  [So,.{r.r3).  5„(„,f,)  •  rl®'(F4)]*  (75) 
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where  n,  i/,  t,  and  k  represent  either  the  horizontal  or  vertical  polarization.  Further¬ 
more,  fi  and  u  are  associated  with  the  receiver  whereas  r  and  k  are  associated  with 
the  transmitter.  By  assuming  gaussian  random  variables  and  considering  only  the 
backscattering  direction,  the  random  term  of  (75)  may  be  written  as  a  sum  of  tw'o 
terms 

<  <?(ri)<?{r2)<?(r3)*(?(f4)*  >  =<  Q(r,)(3{rs)*  ><  <?(r2)(?(r4)*  > 

+  <  Q{rx)Q{uY  ><  Q(rj)(?(r3)’  >  (76) 


Figure  5;  Physical  interpretation  for  Bom  second-order  approximation. 
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Similarly  to  (28),  by  using  the  spectrum  density  of  the  autocorrelation  function  and 
integrating  over  the  lateral  variables,  the  ensemble  average  takes  the  following  form 


<  Ez\r)  >=  d%  j  d.,d^,dz,dz, 

g>(p*ui+«*u)*i  ^  O^p) 

a,7=*.» 

where  the  following  terms  have  been  defined 

=  x,^x„iX*^,x:,,  (78 

Tr;lfa(fce)  =  C^J^(kp)C',W(kp)  (79: 

n!r.'i!;''''(*e)  =  lA(pfci«)  •  &{qku)  a{-rk„)  -  f(sl:i„)l 

ii>(p'ikix.)  ■  i{q%.)  l(-T'ku)  •  ;i(3'*r,„)]*  (80: 

=  (p(p<^ix.)  •  a{qhz)  Q(-rAru)  •  7(3*,^)) 

li>(j>'Jtix.)  •  7(-?'i-ix)  7(r'^ix)  •  K(s'fc,ri)]*  (8I: 

In  the  polar  coordinate  system,  (77)  becomes 

<  EZ\f]  4l\fy  >=  ^Ori^o^i/ dk,  J dzidz2d:idi, 

kp  )xj  >(•*>..+’•*1.1X1  --•V*r„ +»'**.  h» 

L  T:;-!;:o(*>)  ''(^>) / ‘^•3>jdJx3  +  V*^x.)  *(<:..  - 

g-.fl..(x.-.-,)  ^-.3.,(--,-xx)  ^  /^/,fW-’‘'(it,) I  (fAjdJ.j  $(A>.  + 

where  the  integration  over  ^  has  yielded  [3] 

1  /'do  ‘'(A,)  = 

Ti  VO 


(83) 
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.S-i 


Using  the  Cauchy's  theorem,  the  integration  over  3,i  and  J.-j  may  be  carried  out 
The  final  result  for  the  polarimetric  backscatteriug  coefficients  are  obtained  ^3j 


/ 


dk. 


iRes^ikpi  ■rkf.ii' / 


iRes^{k^-k,,0:)  £  Y,  E 


86. 


The  factors  M  and  .V  are  defined  in  (3].  It  is  noted  that  the  integral  ovtr  c.  cannct 
be  solved  aualttically  but  only  numerically.  From  the  integration  over  ;ae  dzimuiaa. 
variable  p  in  (83)  and  (84),  it  is  found  that  =  OkKkv  -  0.  However,  tiie 
depolarization  term  o/,,  is  different  from  zero,  unlike  the  result  obtained  with  the 
isotropic  Bom  first-order  approximation. 
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Results 

a)  Born  first-order  approximation 

With  the  Bom  first-order  approximation,  the  polarimetric  backscattering  coeffi¬ 
cients  have  been  obtained  for  the  two-layer  isotropic  and  anisotropic  random  medium 
configurations.  For  the  anisotropic  random  medium,  the  depolarization  term  a*,  0 

exists  even  with  Born  first-order  approximation.  In  general,  all  the  backscattering 
coefficients  o**,  o*,,  a„.  a**,,,  <Ta„„  are  different  from  zero  for  the  anisotropic 

case.  However,  when  =  V’/  =  0*  but  ei  ^  fi,,  no  depolarization  term  is  obtained  for 
this  uniaxial  untilted  case  and  =  Okv-n  =  0.  Also,  for  V"  =  0  but  V’/  ^  0. 

the  non-zero  depolarization  terms  OkAAr  and  are  obtained.  Finally,  when  V  r  0, 
all  the  backscattering  coefficients  are  different  from  zero.  In  a  similar  manner,  when 
0  =  90°.  implying  that  the  plane  of  incidence  contains  the  scatterers,  there  is  no  de¬ 
polarization  factor,  and  only  three  backscattering  coefficients  are  different  from  zero 

(^Ac  ~  ^AAAr  —  ”  0). 

Fig.  6  illustrates  dependence  of  the  backscattering  coefficients  on  the  angle 
of  incidence  for  <7**,  <r„  and  for  a  two-layer  configuration,  and  these  results 
are  matched  with  experimental  data  for  sea  ice  taken  in  1984  in  Point  Barrow  in 
Alaska  by  [lOj.  The  input  parameters  fj,  ei^,  6,  and  6,  are  obtained  using  the  strong 
fluctuation  theory'  developed  in  (ll].  The  tilt  angle  is  chosen  to  be  =  v/  =  20°, 
and  the  azimuthal  angle  is  assumed  to  be  ^  =  20*.  The  good  correspondence  is 
noted  between  the  experimental  results  and  the  theoretical  predictions.  Fixing  the 
incident  angle  at  ^o.  =  40°  in  Fig.  6,  the  covariance  matrix  for  the  backscattering 
coefficient  (14)  is  shown  in  Fig.  7.  As  can  be  seen,  this  matrix  does  not  contain  any 
zero  elements. 

For  an  earth  terrain  medium  sucii  as  trees  or  grass,  leaves  grow  randomly  in  the 
azimuthal  direction.  This  can  be  modeled  by  assuming  that  one  patch  illuminated 
by  the  radar  is  oriented  in  a  certain  direction  whereas  the  next  patch  is  oriented  in  a 
different  direction  and  so  on.  To  obtain  the  backscattering  coefficients,  the  average 
over  many  patches  should  be  taken.  This  model  is  used  in  the  theoty  by  assuming 
a  fixed  direction  of  the  scatterers  but  modiMng  the  position  of  the  radar  o\Tr  the 
azimuthal  direction.  This  is  shown  in  Fig.  8  for  grass  and  in  Fig.  9  for  trees  where 
both  these  vegetation  fields  are  simulated  by  averaging  the  backscattering  coefficients 
over  o.  The  step  of  averaging  is  Ao  =  1°.  Though  the  tilt  angle  i’  is  different  from 
zero,  there  is  a  depolarization  factor  (oa*  0).  but  OkAAr  =  t^Arrr  =  0.  This  is  due  to 
sjmmetrical  reasons  because  the  contribution  of  these  two  backscattering  coefficients 
at  d  is  cancelled  by  the  one  at  180*  -  p.  Also  in  Fig.  8  and  in  Fig.  9.  experimental 
polarimetric  radar  data  obtained  by  MIT  Lincoln  Laboratory*  are  displayed  [14].  The 
frequency  of  operation  was  35  GHz,  and  the  angle  of  incidence  was  82°  (depression 
angle  of  8°).  The  radar  illuminated  on  a  vegetation  field  consisting  of  either  grass  or 
trees  at  a  range  of  ^  2  km.  A  good  correspondence  between  the  theoretical  predictions 
and  the  experimental  results  is  obtained. 
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Figure  8:  Covariance  matrix  elements  for  grass  region  using  the  average  over  the 
azimuthal  angle  ^  with  =  1*.  Experimental  measurements  and 
theoretical  calculations  are  shown. 
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azimuthal  angle  o  with  So  =  1*.  Experimental  measurements  and 
theoretical  calculations  are  shown. 
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b)  Born  second^order  approximation 

Unlike  the  isotropic  Bora  first-order  approximation,  a  depolarization  factor 
cTit,  ^  0  is  obtained  with  the  isotropic  Born  second-order  approximation  (9).  Fur¬ 
thermore,  it  was  shown  that  o***,  =  =  0.  Hence,  the  covariance  matrix  can  be 

written  from  (11) 

C  =  0  Ok.  0  (86) 

9  <^r.  / 

Although  a  mathematical  proof  has  been  used  to  show  that  five  elements  of  the  covari¬ 
ance  matrix  are  zero,  this  fact  may  also  be  illustrated  from  a  physical  point  of  view. 
Let  us  concentrate  on  the  term  <7„k.  which  is  proportional  to  <  W  •  HV*  >  from 
(11).  By  assuming  that  a  vertically  polarized  electric  field  is  transmitted,  symbolized 
by  V  in  Fig.  10.  the  product  of  the  waves  polarized  in  the  horizontal  and  vertical  di¬ 
rections  at  the  receiver  is  taken.  The  v>zye  hits  a  first  scatterer  at  position  ri  and  then 
hits  a  second  scatterer  at  position  r'^.  At  that  location,  the  electric  field  ^  induces  a 
dipole  which  can  be  decomposed  into  horizontal  and  vertical  components.  From  that 
point,  the  Vr-ave  is  backscattered  and  the  receiver  collects  either  HV*  or  VV.  Due 
to  symmetry  reasons,  the  second  scatterer  may  be  located  at  position  fj.  The  corre- 
spondmg  electric  field  "E"  is  decomposed  into  a  vertical  component  and  a  horizontal 
component  which  is  in  the  opposite  direction  from  the  component  created  by  fj.  At 
the  receiver,  either  HV’  or  VV  is  collected.  In  order  to  obtain  a,,k,,  the  product 
of  these  two  quantities  must  be  taken.  Howver,  the  product  created  by  the  location 
at  jj  wll  have  an  opposite  sign  from  the  one  created  at  the  location  fj.  By  taking 
the  average,  these  two  contributions  cancel  and  o,,*.  =  0  is  obtained.  Similarly  for 
every  point,  there  is  a  symmetric  scatterer  which  will  cancel  its  contribution. 

The  depolarization  factor  is  proportional  to  <  HV  •  HV*  >,  i.e.,  a  vertically 
polarized  electric  field  is  transnutted  and  the  horizontal  component  of  the  scattered 
field  is  received.  Unlike  the  term  <  VV  ■  HV*  >,  both  quantities  HV  and  HV*  due 
to  the  dipole  at  location  F'j  are  in  the  opposite  direction  from  the  ones  created  by  the 
scatterer  at  T'j.  This  implies  that  the  product  of  the  two  quantities,  w'hich  have  the 
same  sign  independent  of  the  position  of  the  second  scatterer  f,  or  rj,  renders  the 
average  backscattering  coefficient  different  from  zero  as  shown  in  Fig.  11. 
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Figure  11;  Interpretation  of  the  term  <  HV  ■  HV‘  >  with  the  Bom  second  order 
approximation. 
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c)  Strong  permittivity  fluctuation  theory 

The  \-ariance  of  permittivity  fluctuations  is  usually  large  in  geophysical  media 
consisting  of  mixtures  of  constituents  with  ver^'  different  properties.  With  the  di¬ 
agrammatic  approach  [l],  the  propagation  and  scattering  in  random  media  can  be 
studied.  The  diagrammatic  approach  leads  to  the  Dyson's  equation  for  the  mean 
field.  It  should  be  noted  that  in  soKnng  Dyson's  equation  with  the  bilocal  approxi¬ 
mation,  both  the  obseiAation  point  and  the  source  point  wthin  the  random  medium 
can  coincide  wth  each  other  in  the  domain  of  integration.  To  take  care  of  the  singu¬ 
lar  nature  of  the  dyadic  Green’s  functions,  a  strong  permittivity  fluctuation  theory 
that  applies  to  both  small  and  large  variances  of  the  permittivity  functions  has  been 
developed  (11). 

The  mean  field  and  the  mean  Green's  function  are  calculated  bj"  replacing  the 
medium  with  an  equivalent  medium  having  an  effective  pennittirity  which  accounts 
for  losses  due  to  both  absorption  and  scattering. 

It  should  be  emphasized  that  the  effect  of  the  random  medium  on  electromag¬ 
netic  waves  is  better  described  with  the  strong  permittivity  fluctuation  theory  than  in 
the  case  of  the  traditional  “weak  permittivity  fluctuations”.  For  the  case  of  a  random 
medium  described  by  a  mixture  of  two  components,  only  the  fractional  volume,  the 
correlation  lengths,  the  permittivity  of  the  two  components,  and  the  frequency  are 
necessary  to  describe  the  effect  of  the  earth  terrain.  The  other  parameters,  that  is 
the  variances  and  the  mean  permittivities,  needed  for  the  Bom  approximation  are 
computed  with  the  strong  permittivity  fluctuation  theory’. 

In  order  to  illustrate  the  results  of  the  strong  permittivity  fluctuation  theory, 
we  compare  <7**  in  Fig.  12  obtained  with  the  weak  permittivity  fluctuation  theory 
compared  with  the  strong  permittivity  fluctuation  theory.  We  use  a  correlation  func¬ 
tion  e.xponentially  decaying  in  all  directions  with  —  U  =  0.3  mm.  The  physical 
parameters  chosen  in  Fig.  12  characterized  a  layer  of  dry’  snow  |13].  In  Fig.  12. 
we  plot  OfcA  in  the  function  of  the  frequency  for  $o,  -  60®.  As  shown,  the  results 
obtained  with  the  strong  permittivity  fluctuation  theory’  are  smaller  than  those  ob¬ 
tained  with  the  weak  permittivity  fluctuation  theon’.  This  difference  gets  larger  with 
increasing  frequencies.  This  is  due  to  the  fact  that  the  losses  due  to  scattering,  which 
increase  with  increasing  frequencies,  are  taken  into  account  in  the  strong  permittivity 
fluctuation  theory. 


Figure  12:  Backscattering  coefficient  <7**  in  function  of  the  frequency  computed 
w,-ith  the  weak  and  strong  pennitthity  fluctuation  theories. 
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